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VISCOUS  FLUID  MOTION  IN  A  SPINNING  AND  NUTATING  CYLINDER 
1.  Accomplishments 

The  working  period  for  this  contract  was  originally  83/08/02  -  85/02/02,  but  has  been 
extended  until  85/04/30  in  order  to  provide  for  completing  the  experimental  setup  for  flow 
visualization.  During  this  working  period,  the  following  personnel  has  been  partly  supported 
under  contract  DAAKl  I-83-K-0011; 

Thorwald  Herbert,  Professor,  Principal  Investigator 

J.  Wallace  Grant,  Assistant  Professor 

Relja  Zivojnovic,  Graduate  Student  (M.  S.  level) 

German  Santos  ,  Graduate  Student  (Ph.D.  level) 

Charlotte  R.  Hawley,  Research  Specialist 

David  Pierpond,  Undergraduate  Student,  has  been  involved  in  the  experimental  work  at  no 
cost.  His  senior  project  is  close  to  completion.  Parts  of  the  theoretical  work  have  been  sup¬ 
ported  by  the  Army  Research  Office  under  Contract  D  AAG29-82-K-0129. 

1.1,  Presentations 

New  research  findings  have  been  reported  at  the  following  conferences: 

(1)  “The  Flow  of  Highly  Viscous  Fluid  in  a  Spinning  and  Nutating  Cylinder,”  1083  Scientific 
Conference  on  Chemical  Defence  Research,  November  14-18,  1083,  Aberdeen  Proving 
Ground,  Maryland. 

(2)  “Highly  Viscous  Fluid  Flow  in  a  Spinning  and  Nutating  Cylinder,”  Second  Army  Confer¬ 
ence  on  Applied  Mathematics  and  Computing,  May  22-25,  1084,  Troy,  New  York. 

(3)  “Instability  of  the  Viscous  Flow  in  aSpinning  and  Nutating  Cylinder,”  ARO  Workshop 
on  Liquid-Filled  Shells,  September  20-21,  1084,  Aberdeen  Proving  Ground,  Maryland. 

(4)  “Instability  of  the  Viscous  Flew-  in  a  Spinning  and  Nutating  Cylinder,”  1084  Scientific 
Conference  on  Chemical  Defense  Research,  November  13-16,  1984,  Aberdeen  Proving 
Ground,  Maryland. 

(5)  An  abstract  of  a  paper  entitled  “On  the  Fluid  Motion  in  Liquid-Filled  Shells,”  has  been 
submitted  for  presentation  at  the  1985  Scientific  Conference  on  Chemical  Defense 
Research,  November  19-22,  1985,  Aberdeen  Proving  Ground,  Maryland. 
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1.2.  Publications 

A  selection  of  results  has  been  published  in  the  following  papers: 

(1)  “The  Flow  of  Highly  Viscous  Fluid  in  a  Spinning  and  Nutating  Cylinder,”  Proc.  1083 
Scientific  Conference  on  Chemical  Defense  Research,  Aberdeen  Proving  Ground,  Mary¬ 
land,  (Eds.  R.  L.  Dimmick,  Jr.  &  M.  Rausa),  Report  CRDC-SP-8-1014,  (1984). 

(2)  “Highly  Viscous  Fluid  Flow  in  a  Spinning  and  Nutating  Cylinder,”  Trans.  Second  Army 
Conference  on  Applied  Mathematics  and  Computing,  Troy,  New  York,  ARO  Report  85-1 
(1985). 

(3)  “On  the  Viscous  Roll  Moment  in  a  Spinning  and  Nutating  Cylinder,”  Proc.  1984 
Scientific  Conference  on  Chemical  Defense  Research,  Aberdeen  Proving  Ground,  Mary¬ 
land,  to  be  published. 

(4)  “Viscous  Fluid  Motion  in  a  Spinning  and  Nutating  Cylinder,”  submitted  for  publication 
in  Journal  of  Fluid  Mechanics. 

Copies  of  the  papers  (1),  (2)  and  of  the  manuscripts  (3)  and  (4)  are  attached  as  Appen¬ 
dices  A.l  -  A.4. 

2.  Technical  Discussion 

2.1.  Introduction 

It  is  well-known  that  spin-stabilized  shells  carrying  liquid  payloads  can  suffer  dynamical 
instability.  For  cylindrical  cavities  and  low  viscosity  of  the  liquid,  the  instability  due  to  basically 
inviscid  inertial  waves  can  be  predicted  by  the  Stcwartson-Wedemeyer  theory  (Stewartson  1959; 
Wedemcyer  1966).  This  theory  rests  on  the  boundary-layer  approach  and  is,  therefore,  res¬ 
tricted  to  the  range  of  sufficiently  large  Reynolds  numbers.  The  instability  of  certain  shells  like 
the  XM761  (D’Amico  1977;  1978),  however,  escapes  such  a  prediction  and  is  also  dis¬ 
tinguished  in  character  owing  to  the  rapid  loss  in  spin  rate.  Experiments  with  a  full-scale 
liquid-filled  cylinder  (Miller  1982)  and  subsequent,  field  tests  (D’Amico  &  Miller  1979)  estab¬ 
lish  that  this  new  flight  instability  is  most  pronounced  for  liquid  fills  of  very  high  viscosity. 

We  have  conducted  an  analysis  of  this  problem  in  order  to  support  the  ongoing  experi¬ 
ments  and  to  independently  obtain  insight  into  the  anatomy  of  the  flow  phenomena.  The  initial 
steps  of  this  analysis  are  reported  elsewhere  (Herbert  1982):  evaluation  of  the  experimental 
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data  base,  dimensional  analysis,  scaling  aspects,  governing  equations,  and  discussion  of  various 
simplifying  assumptions,  Two  observations  in  this  earlier  work  led  to  the  approach  discussed  in 
the  following.  First,  if  the  despin  (negative  roll)  moments  (Miller  1982)  and  void  observations 
(Miller  1981)  are  correlated  with  the  Reynolds  number  Re,  at  least  three  regions  can  be  dis¬ 
tinguished.  At  low  Re,  the  despin  moment  increases  proportional  to  Re,  and  the  void  in  an 
incompletely  filled  cylinder  is  parallel  to  the  spin  axis.  This  suggests  a  simple  fluid  motion  that 
is  essentially  independent  of  the  axial  coordinate,  except  in  the  neighborhood  of  the  end  walls. 
In  a  middle  range  of  Re,  the  despin  moment  assumes  a  maximum,  and  a  wavy  distortion  of  the 
void  seems  to  indicate  a  cellular  structure  of  the  fluid  motion.  This  cellular  motion  can,  in 
principle,  originate  from  hydrodynamic  instability  of  the  basic  flow  with  respect  to  axially 
periodic  disturbances.  At  still  higher  Reynolds  numbers,  the  despin  moment  decreases  with 
increasing  Re  in  a  manner  not  clearly  defined  by  the  few  available  data  points.  The  void  obser¬ 
vations  indicate,  however,  that  the  motion  ultimately  becomes  turbulent. 

Tile  second  observation  is  the  appearance  of  the  nutation  rate  and  angle  as  a  small  param¬ 
eter  in  the  equations  for  the  deviation  from  solid-body  rotation.  The  forcing  term  due  to  nutA- 
lion  can  be  considered  small  enough  for  linearization  of  the  equations  in  the  situations  of  prac¬ 
tical  interest. 

Consequently,  our  research  focused  on  three  topics.  First,  theoretical  analysis  of  "simple” 
fluid  motions  at  low  Reynolds  number  that  satisfy  the  linearized  equations  for  the  deviation 
from  rigid  body  rotation.  Second,  analysis  of  inviscid  and  viscous  inertial  modes  and  their 
relevance  for  the  occurrence  of  a  maximum  despin  moment  at  moderate  Reynolds  numbers. 
Third,  the  design  of  a  small-scale,  low-cost  experiment  for  visualization  of  the  interior  fluid 
motion.  The  results  of  these  efforts  are  discussed  in  sections  2.2,  2.3  and  2.4,  respectively. 

2.2.  The  Deviation  from  Solid  Body  Rotation 

A  formal  analysis  of  the  equations!  for  the  deviation  of  the  velocity  field  from  solid  body 
rotation  suggests  that  (1)  the  equations  can  be  linearized  without- introducing  major  errors,  (2) 
at  low  Reynolds  numbers  Rc  ,  the  velocity  held  is  independent  of  the  axial  direction  over  a 
considerable  part  of  the  relatively  long  ( q  =  1.3)  cylinder,  and  (3)  the  (low  is  in  the  axial 
direction  and  turns  at  die  ends  Application  of  these  eonccptiiJ  assumptions  turned  out  very 


t  Detailed  equations  arc  given  in  Appendix  A  1 
•  We  use  tli e  notation  introduced  in  Appendix  A  .1 
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fruitful.  A  detailed  description  of  the  results  and  comparison  with  computational  and  experi¬ 
mental  data  is  given  in  the  paper  '“Viscous  Fluid  Motion  in  a  Spinning  and  Nutating  Cylinder” 
(Appendix  A  .4),  submitted  for  publication  in  the  Journal  of  Fluid  Mechanics.  A  sample  pro¬ 
gram  with  the  relevant  subroutines  for  reproducing  the  theoretical  results  is  given  in  Appendix 
B.  Here,  we  report  only  the  main  conclusions. 

The  model  ol  a  two-dimensional  unidirectional  flow  in  a  finite  segment  of  an  infinite 
cylinder  yields  the  solution  of  the  linearized  equations  in  analytical  form.  The  disregard  of  the 
end  walls  has  some  obvious  consequences:  the  turning  flow  near  the  ends  and  the  associated 
contributions  of  pressure  and  shear  stresses  to  the  moments  cannot  be  obtained  from  this 
model.  Nevertheless,  we  gather  understanding  as  well  as  quantitative  information.  The  velocity 
field  of  the  core  flow  agrees  well  with  computational  results  (Vaughn  et  al.  1985)  for  low  Rey¬ 
nolds  numbers.  The  analytical  result  is  an  evident  example  for  the  formation  of  boundary 
layers.  The  cc  e  flow  can  be  utilized  as  a  basic  flow  in  studies  of  hydrodynamic  instability  with 
respect  to  cellular  motions.  The  parametric  excitation  of  such  cells  by  the  azimuthally  periodic 
deviation  has  been  discussed  by  Herbert  (  1984).  The  core  flow  also  represents  the  lowesteorder 
approximation  to  the  solution  of  the  nonlinear  equations  and  can  be  extended  by  higher-order 
terms. 

The  roll  moment  agrees  well  with  measured  and  computed  values,  and  can  also  be  found 
at  Reynolds  numbers  too  large  for  successful  numerical  simulations.  The  roll  moment  ori¬ 
ginates  from  Coriolis  forces.  While  the  direct  calculation  of  the  yaw  moment  suffers  from 
neglecting  the  pressure  contribution,  the  yaw  moment  can  be  found  from  the  roll  moment 
using  the  relations  given  by  Murphy  (198  1,  1985).  The  pitch  moment  remains  an  open  issue. 
The  average  rate  of  change  of  temperature  is  found  to  be  proportional  to  roll  moment  and  spin 
rate.  This  estimate  needs  further  verification  once  more  detailed  data  become  available. 

The  simple  form  and  scaling  relations  of  our  results  provide  guidance  for  sorting  and 
evaluating  the  experimental  data  base.  The  results  also  suggest  various  improvements  in  the 
experimental  procedures.  First,  the  changes  in  temperature  and  viscosity  should  be  carefully 
monitored  With  the  effective  viscosity  known,  a  closer  agreement  between  theory  and  obser¬ 
vation  is  to  be  expected.  Second,  the  yet  neglected  variation  of  the  roll  moment  with  the  spin 
rate  is  considered  relevant  and  in  fact  provides  the  roll  moment  in  some  range  of  ueynolds 
numbers  Instead  of  recording  the  roll  moment  as  a  function  of  Re  by  using  numerous  viscosi¬ 
ties  at  fixed  spin  rate,  very  similar  data  can  be  generated  by  varying  the  spin  rate  for  a  few 
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fluids.  For  directing  the  research  efforts  within  this  project,  it  has  been  most  revealing  that  the 
characteristic  variation  of  roll  moment  versus  Reynolds  number,  in  particular  the  sharp  max¬ 
imum  at  Re  »  10  is  a  property  of  the  unidirectional  model  flow.  This  result  contradicted  the 
earlier  working  hypothesis  which  attributed  the  occurrence  of  this  maximum  to  hydrodynamic 
instability  and  the  onset  of  cellular  motions. 

2.3.  Cellular  Structure  of  the  Motion 

Although  not  as  relevant  to  the  moments  as  earlier  thought,  the  cellular  motion  at 
surprisingly  low  Reynolds  numbers  is  in  itself  an  interesting  physical  phenomenon,  in  our 
analysis,  we  superpose  to  the  steady  basic  flow  v,  =  (O.r.ti,),  p ,  disturbances  v  =  (u,v,w),  p 
sufficiently  small  for  linearization.  Substitution  into  the  governing  equations  and  neglect  of 
products  between  disturbances  and  terms  of  order  0(e2)  provide  the  following  stability  equa¬ 
tions: 


dv  dv, 
dt  dtp 


2(1  +  t,)v  + 

dr 


Re 


ti  2  dv  | 

r2  r2  c ><j>  ■ 
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Three  groups  of  terms  have  been  separated  by  square  brackets  in  tiiese  ■  ■<  The 

first  group,  if  set  to  zero,  represents  the  equations  for  in  viscid  inertia,  inodes 
~-cxp(im<P  -+  iq:  -i-.7<),  where  m  is  the  (integer)  azimuthal  wavenumber,  a  the  axial 
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wavenumber,  and  er  =  u ,  +  iff,  provides  tlie  amplification  rate  a  r  and  frequency  <r,.  Usually, 
an  equation  for  the  pressure  is  used  for  obtaining  the  analytical  solution  (Stewartson  1959). 
We  have  derived  an  alternative  system  in  terms  of  u  ,v  and  applied  a  spectral  collocation 
method  to  be  used  for  more  general  cases  in  order  to  check  the  numerical  results  against  the 
exact  values. 

The  programs  are  designed  to  provide  results  for  arbitrary  values  of  the  azimuthal  and 
axial  wavenumber.  A  typical  spectrum  for  »i  =  1 ,  a  =  1  is  shown  in  figure  1.  In  the  complex 
s-plane,  i  =  a  +im,  stable  (s,  <  0)  eigenvalues  are  located  to  the  left  of  the  vertical  «,-axis, 
while  eigenvalues  in  the  right  half-plane  indicate  instability.  In  the  inviscid  case,  the  eigen¬ 
values  are  located  on  the  s,-axis  (a,  =  0)  and  are  neutral.  The  eigenfunctions  associated  with 
the  eigenvalue  of  maximum  frequency  are  shown  in  figure  2.  Note  that  only  the  no- 
penetration  condition  for  the  radial  velocity  component  is  satisfied,  while  v  and  tv  slip  at  the 
cylinder  walls. 

The  second  group  of  terms  in  the  stability  equations  is  multiplied  by  1  /Re  and  represents 
the  viscous  correction  to  the  inertial  modes.  By  eliminating  w  and  p,  a  system  of  ordinary 
differential  equations  for  u  and  e  has  been  derived.  Due  to  the  higher  order  of  this  system,  all 
boundary  conditions  can  be  satisfied.  Programs  have  been  developed  for  calculating  spectra  of 
(complex)  eigenvalues,  for  tracing  single  eigenvalues  as  a  function  of  Re,  or,  and  for  obtaining 
the  eigenfunctions.  At  high  Re,  the  results  follow  the  trends  predicted  by  asymptotic  theories. 
Our  analysis,  however,  also  covers  the  range  of  low  Reynolds  numbers,  where  the  inertial 
modes  suffer  rapid  decay  (»P<0). 
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Spectra  of  eigenvalues  for  viscous  inertial  inodes  at  m  =  1 ,  o  =  1  in  an  infinite  cylinder 
are  shown  in  figure  3  for  Re  =  1000  and  in  figure  1  for  Re  =  100  All  eigenvalues  move  to 
the  left  as  Re  decreases,  indicating  stabilization  of  the  modes.  The  eigenfunctions  for  the  least 
stable  mode  at  Re  =  100  arc  shown  in  figure  5.  Note  that  all  velocity  components  vanish  at 
the  cylinder  wall. 
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'flie  most  interesting  aspect  of  the  stability  equations  is  the  third  group  of  terms.  The 
coefficients  in  this  group,  v.,  t,,  and  t  are  of  order  0(  i)  and  periodic  in  <?.  The  periodicity  in 
r p  leads  to  a  coupling  of  the  mode  equations  for  m  and  in  ±  1,  and  may  cause  primary  resonance 
between  inertial  modes.  In  view  of  viscous  damping,  this  resonance  is  likely  to  occur  as  t 
exceeds  a  critical  value  that  decreases  as  Rc  increases. 
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The  analysis  of  this  parametrical  instability  has  been  prepared  by  solv>ng  the  full  set  of 
stability  equations  with  the  coupling  terms  between  modes  artificially  set  to  2ero.  A  spectrum 
of  eigenvalues  in  the  <r-plane  for  tn  =  0,  a  =  1  and  Re  —  100  is  shown  in  figure  6.  This 
spectrum  consists  of  the  superposition  of  the  spectra  for  m  —  0,  ±  1  of  the  viscous  equations. 
As  the  coupling  terms  are  taken  into  account,  the  spectrum  drastically  changes,  as  shown  in 
figure  7.  Most  remarkable  is  the  appearance  of  a  small  number  of  eigenvalues  in  the  unstable 
right  half-plane.  These  eigenvalues  are  either  real  or  appear  as  complex  conjugate  pairs.  This 
reflects  the  fact  that  the  physical  solution  must  be  real.  Instability  at  lower  values  of  the  axial 
wavenumber  a  can  be  found  at  low  Reynolds  number.  Figure  8  shows  an  example  of  instabil¬ 
ity  at  Re  =  20.  A  single  pair  of  complex  conjugate  eigenvalues  appears  in  the  unstable 
domain;  the  amplification  rate  decreases  with  the  Reynolds  number. 

The  analysis  of  the  parametric  instability  of  the  azimuthally  periodic  flow  is  rather  costly 
in  terms  of  computer  time.  Within  the  framework  of  this  contract,  we  were  not  able  to  com¬ 
plete  a  systematic  parameter  study  that  could  provide  a  lower  limit  or  critical  Reynolds  number 
for  the  onset  of  cellular  motion;  neither  were  experimental  data  available  for  a  comparison  of 
eigenfunctions.  It  has  been  shown,  however,  that  cellular  motions  can  appear  at  low  Reynolds 
numbers  as  a  result  of  resonant  coupling  beetween  inertial  modes. 
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2.4.  Flow  Visualization 

Although  theoretical  (Herbert,  Appendix  A  A)  and  computational  (Vaughn  et  al.  1985) 
work  provides  some  insight  into  the  interior  fluid  motion,  the  nature  of  the  phenomena 
remains  largely  in  the  dark.  This  is  especially  true  for  the  range  of  medium  and  high  Reynolds 
numbers  where  finite-amplitude  cellular  motions  and,  ultimately,  turbulence  are  expected  to 
occur.  Tli is  range  is  barely  within  the  scope  of  computational  methods  nor  can  it  be  fully 
explored  with  the  theoretical  means  of  sections  2.2,  2.3. 

Previous  experiments  (Miller  1981)  using  a  partially  filled  full-scale  cylinder  revealed 
some  axial  nonuniformitv  of  the  flow  at  higher  Re  without  showing  details  of  the  flow  field. 
Later  attempts  to  use  flow  tracers  (Miller  .t  Oberkampf,  personal  communication)  had  little 
success  due  to  the  high  spin  rates  (accelerations)  combined  wi»h  minute  density  dilTerenccs 
between  working  fluid  and  tracer  particles.  Even  carefully  centrifuged  and  selected  particles 
failed  to  follow  the  liquid  path,  probably  due  to  changes  of  i<  mperature  during  the  run. 
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Attempt*  to  employ  laser-induced  fluorescence  (Miller  1084)  were  partly  successful  after  chang¬ 
ing  the  time  scale,  i.e,  to  lower  spinrate,  nutation  rate,  and  viscosity  at  fixed  values  of  the 
dimensionless  parameters.  These  efTorts  have  been  discontinued,  however,  due  to  continuing 
lighting  problems. 

Evaluation  of  the  experimental  attempts  to  visualise  the  fluid  flow  clearly  reveals  the 
extreme  full-scale  conditions  as  evil.  Conclusive  experiments  can  be  conducted  by  exploiting 
the  principals  of  dynamical  similarity  and  appropriate  scaling  laws.  These  aspects  have  been  dis¬ 
cussed  in  earlier  work  (Herbert  1082)  and  extended  by  the  analysis  in  Appendix  A. 4. 

Between  the  three  reference  quantities,  radius  a,  spinrate  w,  and  density  p  for  length, 
time,  and  mass,  respectively,  the  density  of  dilTerent  fluids  offers  little  variability.  However, 
length  and  time  scale  can  be  easily  changed.  For  dynamical  similarity,  the  following  dimension¬ 
less  quantities  must  be  fixed: 


\  —  c  /a 


aspect  ratio 
nutation  angle 
frequency 


r  =  n  /w  frequency 

Re  =  pwa2/^  Reynolds  number 


The  nutation  angle  must  remain  the  same  in  a  scaled  setup.  Radius  a  and  half-length  c  of  the 
cylinder  must  be  scaled  by  the  same  factor  in  order  to  keep  the  aspect  ratio  fixed.  A  second 
factor  can  be  applied  to  both  spinrate  w  and  nutation  rate  0  ,  in  order  to  preserve  the  fre¬ 
quency.  Keeping  Re  fixed  requires  changing  the  kinematic  viscosity  v  =-  /j ,/p  by  the  same  fac¬ 
tor  as  u I  a'1.  Since  the  desired  tendency  is  toward  smaller  radii  and  spinrates,  we  require  less 
viscous  fluids  than  those  used  in  the  full-scale  experiments.  Such  fluids  arc  easy  to  find. 

It  is  obvious  that  the  main  thrust  of  an  experiment  may  require  specific  optimum  condi¬ 
tions.  Flow  visualization  requires  low  velocities,  i.e.  low  values  of  wa.  Measurements  of 
moments  require  optimum  values  of  wsaJ.  Minimizing  the  rate  of  change  of  temperature 
requires  a  minimum  of  w3u2.  A  good  setup  for  flow  visualization,  therefore,  may  produce 
moments  in  a  hardly  measurable  range. 

The  goal  of  our  efTorts  was  to  show  that  a  low-cost  device  (s*$0UU)  can  be  designed  for 
flow  visualization.  Details  had  to  be  kept  simple.  Accuracy  and  convenience  had  to  comprom- 
inc.  The  results  of  these  efforts  arc  t>hown  in  Figures  0  and  10. 
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A  one-inch  inner  diameter  cylinder  of  aspect  ratio  X  =  4.3  is  used.  The  cylinder  is  cut 
from  a  pyrex  glass  tube  with  the  inner  diameter  accurate  within  1/5000  inch,  but  with  varying 
wall  thickness  that  affects  the  optical  quality.  The  cylinder  is  filled  with  mixtures  of  water  and 
glycerin.  The  mixing  ratio  is  used  to  vary  viscosity.  On  top,  the  cylinder  is  closed  with  a 
screwed-in  plastic  plug.  A  center  hole  allows  access  to  the  interior,  especially  for  removing  air 
bubbles,  The  hole  can  be  closed  using  a  toothpick. 

The  cylinder  is  glued  to  a  drive  plug  and  axis  machined  from  a  single  piece  of  aluminum. 
The  one-sided  support  allows  easy  (optical)  access  to  the  cylinder  and  permits  using  cylinders  of 
different  length.  One-sided  support  is  affordable  due  to  the  moments  being  approximately  five 
orders  of  magnitude  smaller  than  in  the  full-scale  experiments.  The  axis  is  twice  supported  by 
ball  bearings.  The  cylinder  and  shaft  are  driven  via  timing  belts  over  exchangeable  sets  of  pul¬ 
leys.  The  most  expensive  piece  of  the  spin  arrangement  is  a  <2lV  d.c.  motor  with  sufficient 
torque  in  the  range  of  500  -  5000  rpm.  Motor  and  cylinder  support  are  mounted  to  an  alumi¬ 
num  frame  that  can  be  inclined  to  the  vertical  axis  by  5,  10,  15  and  20 ‘using  different  support 
holes  and  struts. 

The  horizontal  support  plate  is  machined  to  leave  the  center  position  free  for  access  and  is 
screwed  to  a  commercial  record  player  (Garrard  model  776).  The  plate  can  be  offset  in  order  to 
align  the  center  of  mass  of  the  cylinder  with  the  nutation  axis.  The  record  player  provides 
nutation  rates  of  33,  45,  and  78  rpm.  The  hollow  axis  is  only  utilized  to  provide  power  to  the 
spin  motor.  A  10-contact  slip  ring  of  high  quality  has  been  kindly  donated  by  Poly-Scicntific, 
Blacksburg,  VA,  but  has  been  saved  for  more  sophisticated  experiments.  A  nail  with  a  smooth 
top  and  a  brush  fixed  to  the  turntable  proved  sufficient  for  transmitting  a  single  voltage. 

The  remaining  components  of  the  experiment  are:  a  Heathkit  regulated  power  supply  for 
the  spin  motor,  astrobelight  for  controlled  pulsed  lighting,  and  suitable  flow  tracers. 

As  flow  tracers  we  use  Afflair  100  Silver  Pearl,  kindly  donated  by  EM  Chemicals, 
Hawthorne,  NY.  The  material  consists  of  very  fine  and  shiny  plastic  platelets  commercially 
used  for  cosmetic  purposes.  Although  of  specific  weight  different  from  that  of  the  fluid,  the 
low  accelerations  in  the  scale  model  permit  practically  buoyant  behavior  of  the  platelets  over 
considerable  time.  The  strobcliglil  ( General  Radio  Strobotac)  with  adjustable  frequency  is  used 
for  lighting  as  well  as  for  measuring  the  spin  rate  of  the  cylinder. 

At  the  slow  time  scale  of  the  experiment,  the  fluid  motion  can  be  visually  inspected  while 
running  the  apparatus.  At  high  viscosities,  the  apparatus  can  also  be  suddenly  stopped,  with  the 
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flow  tracers  “frozen”  in  the  resting  fluid.  The  platelets  align  with  surfaces  of  constant  shear. 
Therefore,  by  manually  rotating  the  cylinder' forth  and  back,  the  three-dimensional  structure  of 
the  field  can  be  inspected.  This  crude  observation  is  very  helpful  in  developing  the  visualiza¬ 
tion  technique.  A  detailed  account  of  the  technique  (appropriate  particle  density,  pitfalls  such 
as  the  history  of  particle  distribution  and  alignment)  will  be  given  by  Pierpont  ( 1985) . 

A  typical  view  of  the  “frozen”  pattern  is  shown  in  figure  11.  Note,  however,  that  the 
stereoscopic  view  by  eye  reveals  the  spatial  distribution  of  the  particles  which  is  here  projected 
into  a  single  plane.  The  photograph  also  shows  some  undesirable  reflections  from  the  cylindri¬ 
cal  surfaces.  Visualization  of  the  frozen  pattern  can  be  essentially  improved  by  using  a  light 
sheet  parallel  to  the  spin  axis.  Sheet  lighting  also  enables  photographic  recording  of  the  flow 
structure  while  the  apparatus  is  in  operation. 

A  continuous  light  sheet  is  produced  by  a  Speclia  Physics  model  120  (15  m\V)  helium- 
neon  laser  and  a  cylinder  lens.  In  order  to  avoid  the  need  for  accurately  firing  the  camera  (35 
mm  Pcntax  with  50  mm  lens)  at  a  certain  time,  a  cylindrical  card  board  screen  with  a  vertical 
slot  and  a  00’  offset  opening  is  fixed  to  the  circumference  of  the  turntable.  The  shutter  is 
manually  opened  and  closed  after  the  laser  sheet  of  light  (lashed  through  the  slot.  Some  photo¬ 
graphs  taken  with  the  apparatus  in  motion  are  shown  in  the  figures.  The  figures  show  the  pair 
tern  in  the  x,z-plane  (the  plane  spanned  by  spin  axis  and  nutation  axis)  as  seen  at  <t>  ==  0’  in 
the  y-direction.  The  series  of  figures  12, 13,14  is  for  different  nutation  rates  of  33,  45,  78  rpm, 
respectively  at  otherwise  fixed  parameters.  The  figures  reveal  the  changing  pattern  and  formar 
tion  of  cells  at  a  Reynolds  number  Re  ^  40.  At  slightly  higher  Re,  additional  cells  appear  as 
shown  in  figure  15.  At  the  high  value  of  Re  8000,  the  flow  is  highly  unsteady  and  irregular 
(turbulent)  with  a  superposed  large-scale  cellular  motion.  This  cellular  motion  accumulates  the 
platelets  in  streaks  around  the  cylinder  wall,  as  shown  in  figure  10.  Number  and  position  of  the 
streaks  depends  sensitively  on  the  parameters,  see  figure  17.  The  distortion  of  the  bright  line 
near  the  axis  is  very  similar  to  the  distortion  oT  the  void  in  Miller's  (1081)  photographs  at  low 
viscosity. 

Specification  of  accurate  Reynolds  numbers  suffers  from  some  uncertainty  in  monitoring 
and  measuring  the  wide  range  of  viscosities  for  the  hygroscopic  water-glycerin  mixtures  exposed 
to  uncontrolled  thermal  conditions,  Falling-sphere  viscomctry  is  the  easy  way  out,  but  requires 
a  whole  variety  of  spheres,  different  in  diameter  and  specific  gravity.  To  within  this  uncer¬ 
tainty.  however,  the  figures  clearly  reveal  the  cellular  structure  of  the  flow  and  the  changes  of 
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the  structure  as  Reynolds  number  and  rale  increase.  Perhaps  the  most  striking  result  of  this 
visual  study  of  the  flow  structure  u  the  manifold  of  pattern  at  higher  Reynolds  numbers.  A 
systematic  analysis  of  these  patterns  cannot  be  conducted  within  the  ongoing  senior  project. 
Neither  can  a  coarse  sampling  of  pattern  in  the  y, c-plane  be  completed  before  submission  of 
this  report.  However,  the  feasibility  of  flow  visualization  with  simple  means  by  proper  sealing 
has  been  clearly  demonstrated.  A  manifold  of  cellular  motions  in  the  laminar  range  has  been 
observed  as  well  as  turbulent  flow  at  high  Reynolds  numbers. 

Some  improvements  are  suggested  from  experience  with  the  present  experiment: 

(1)  More  accurate  measurements  of  viscosity  and  spinrate  in  order  to  reproduce  experimental 
conditions. 

(2)  Redesign  of  the  spin  drive  plug:  a  steel  axis  should  be  fixed  to  the  bearings  with  an  easy- 
to-detach  cylinder  bottom  plug  snug-fit  to  the  axis. 

(3)  Refined  laser-light  sheet  (thinner,  plane  can  be  rotated)  and  more  sophisticated  shutter 
release  for  the  camera  (controlled  by  timing  pulse  at  proper  turntable  position). 

(4)  Higher  precision  of  the  cylinder  for  improved  optical  quality. 

(5)  Replacement  of  the  record  player  by  a  turntable  with  continuously  variable  rate  of  rotation 
and  heavier  duty  bearings. 

Various  of  these  improvements  can  be  achieved  with  minor  efforts  and  expenses,  but  not 
within  the  time  frame  of  this  contract.  Moreover,  we  are  confident  that  using  photochromic 
dye  exposed  to  a  pulsed  laser  beam  would  allow  visualization  and  -  with  some  effort  -  quantita¬ 
tive  analysis  of  velocity  profiles 
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THE  FLOW  OF  HIGHLY  VISCOUS  FLUID 
IN  A  SPINNING  AND  NUTATING  CYLINDER 


Thorwald  Herbert 

Department  of  Engineering  Science  and  Mechanics 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  Virginia  24061 


Abstract 

Evaluation  of  the  sparse  experimental  data  and  analysis  of  the  equations 
for  the  fluid  motion  In  a  spinning  and  nutating  cylinder  suggest  a  theoretical 
analysis  In  three  steps:  (1)  analysis  of  the  viscous  flow  In  an  infinitely 
long  cylinder,  (2)  hydrodynamic  stability  analysis  of  this  basic  flow,  and  (3) 
study  of  the  end  effects.  A  status  report  Is  given  on  the  first  of  these 
tasks.  The  purely  axial  basic  flow  Is  governed  by  ordinary  differential 
equations.  Analytical  and  numerical  solutions  show  the  different  character  of 
this  flow  at  low  and  high  Reynolds  numbers.  The  resulting  moments  are 
consistent  with  experimental  data. 


Introduction 

Spin-stabilized  projectiles  with  liquid  payloads  can  experience  a  severe 
flight  Instability  that  is  characterized  by  a  rap  yaw  angle  growth  and  a 
simultaneous  loss  In  spin  rate.  Experiments  with  u  full-scale  liquid-filled 
cylinder  (Miller  1982)  have  shown  that  this  Instability  originates  from  the 
Internal  motion  of  fluids  In  the  range  of  high  viscosities.  We  have  Initiated 
a  theoretical  analysis  of  this  problem  In  order  to  support  the  ongoing  experi¬ 
ments  and  to  obtain  Independently  Insight  Into  the  flow  phenomena.  The 
Initial  steps  of  our  approach  are  reported  elsewhere  (Herbert  1982):  evalua¬ 
tion  of  the  experimental  data  base,  dimensional  analysis,  scaling  aspects, 
governing  equations,  and  discussion  of  various  simplifying  assumptions.  Two 
observations  In  this  earlier  work  led  to  the  results  presented  here.  First, 

If  the  despin  moment  data  and  void  observations  (Miller  1981)  are  correlated 
with  the  Reynolds  number,  at  least  three  regions  can  be  distinguished.  At  low 
Reynolds  numbers  Re,  the  despin  moment  Increases  proportional  to  Re.  and  the 
void  In  an  Incompletely  filled  cylinder  Is  parallel  to  the  spin  axis.  In  a 
middle  range,  a  maximum  of  the  despin  moment  exists  that  seems  to  be  associ¬ 
ated  with  the  onset  of  a  void  distortion  due  to  a  cellular  motion  of  the 
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fluid.  At  higher  Reynolds  numbers,  the  despin  moment  decreases  In  a  manner 
that  is  not  clearly  defined  by  the  few  available  data  points.  The  void 
observations  Indicate,  however,  that  the  motion  ultimately  is  turbulent.  In 
the  theoretical  approach,  the  Initial  Increase  In  the  despin  moment  Is 
attributed  to  a  simple  basic  flow  that  can  be  studied  disregarding  the  end 
effects.  The  occurrence  of  a  cellular  motion  at  higher  Reynolds  numbers  and 
tne  reduced  despin  moment  can.  In  principle,  be  associated  with  a  hydrodynamic 
Instability  of  the  basic  flow.  Accounting  for  the  end  effects  should  provide 
the  final  polish  of  the  results.  The  second  observation  was  the  appearance  of 
a  parameter  in  the  equations  for  the  deviation  from  solid  body  rotation  that 
can  be  considered  small  enough  for  linearization  in  the  situations  of  practi¬ 
cal  and  experimental  Interest. 

In  the  following,  we  describe  the  development  of  a  simple  system  of 
equations  for  the  basic  flow  and  discuss  asymptotic  solutions  end  numerical 
results.  A  comparison  Is  made  with  computer  simulations  of  the  flow  (Vaughn 
1983)  and  with  experimental  data  (Miller  1982). 


The  Basic  Flow 

We  consider  the  motion  of  a  fluid  of  density  p  and  viscosity  u  In  a 
cylinder  of  radius  a  and  length  2c.  The  cylinder  rotates  with  the  spin  rate  u 
about  Its  axis,  the  z-axls.  The  z-axls  Is  Inclined  by  the  nutation  angle  e 
with  respect  to  the  Inertial  2-axis  and  the  2,z-plane  rotates  with  the  nuta¬ 
tion  rate  n  about  the  Z-axls.  The  two  rotation  axes  Intersect  In  the  center 
of  mass  of  the  cylinder.  In  contrast  to  the  experimental  procedures  (Miller 
1982),  we  consider  w»0,  fl>0,  and  o>0  as  constant.  We  describe  the  motion  by 
the  Navler-Stokes  equations,  written  In  the  aeroballlstlc  or  nutating  coordi¬ 
nate  system  x.y,2,  where  x  Is  normal  to  z  In  the  Z,z-plane: 

DV 

+  22  *  X*  4  S  *  (£  *  £)]  *  -vPn  +  n^Yn 

v»V„  -  0. 

~n 

Y_  Is  the  velocity  measured  In  the  nutating  system,  pn  the  pressure,  and  £ 
the  position  vector.  The  body  force  due  to  gravity  Is  neglected.  Equations 
(1)  are  subject  to  the  no-sllp  and  no-penetration  conditions  at  the  cylinder 
walls. 

In  the  next  step,  we  split  the  velocity  and  pressure  fields  according  to 

in  "  is  +  V  Pn  “  Ps  +  Pd  (2) 


where  Ys, p.  describe  the  state  of  pure  solid  body  rotation,  whereas 
describe  the  deviation  from  solid  body  rotation.  It  Is  obvious  that’ )LaS, 
Pd=0  If  either  one  of  the  following  conditions  Is  satisfied:  u*0,  n«0,  e-0 
or  y—  (solid  fill). 

The  equations  for  are  written  in  terms  of  non-dimensional  quantities, 
using  a,  w,  and  p  for  scaling  length,  time  and  mass.  Note  that  this  choice  Is 
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ambiguous  (Herbert  1982)  and  excludes  the  case  of  no  spin.  The  problem  then 
depends  on  four  non-dimensional  parameters: 


\  -  c/e 
a  -  sine 
t  ■  n/w 
Re  ■  pwa2/p 


aspect  ratio 

frequency 
Reynolds  number. 


The  aspect  ratio  Is  only  contained  In  the  boundary  conditions. 


Finally,  cylindrical  coordinates  r,*,z  are  Introduced.  The  equations  for 
the  (non-dimensional)  components  vr,v.,v_  of  the  velocity  deviation  and  the 
pressure  deviation  p*  are  given  as  equations  (4.22)  by  Herbert  (1982).  The 
most  remarkable  fact  about  these  equations  Is  the  appearance  of  a  force  term 
2?0rco$t  in  the  z-momentum  equation.  If  this  term  vanishes,  the  equations 
support  only  a  trivial  solution.  It  is  obvious,  therefore,  that  the  deviation 
velocity  Is  of  order  0(to).  At  close  analysis,  to  ■  (fl/wjslne  turns  out  to 
be  a  rather  small  parameter.  Even  a  conservative  estimate  with  n<500  rpm, 
w>3000  rpm  and  e<20°  provides  values  to<0.057.  It  seems  well  Justified, 
then,  to  linearize  the  equations  In  to.  It  Is  also  worth  noting  that  this 
linearization  Imposes  no  restriction  on  the  Reynolds  number  Re. 


The  resulting  system  of  linear  equations  Is  still  quite  difficult  to 
solve,  mainly  owing  to  the  boundary  conditions.  Use  of  the  boundary  layer 
approximation  would  simplify  the  task,  but  seems  Inappropriate  In  the  Inter¬ 
esting  range  of  Reynolds  numbers.  The  equations  suggest,  however,  that  a 
solution  can  be  found  when  the  boundary  conditions  at  the  end  walls  of  the 
cylinder  are  relaxed.  As  a  first  step,  therefore,  a  steady  flow  (basic  flow) 
Is  sought  disregarding  the  end  walls,  l.e.  for  an  infinitely  long  cylinder. 
Despite  other  Initial  thoughts,  the  force  term  In  the  z-momentum  equation 
causes  a  deviation  velocity  that  Is  purely  axial, 

Vd  ■  (0,  0,  «avz),  pd  •  r).  (3) 

It  Is  also  consistent  with  the  equations  to  assume  a  solution  In  the  form 


where 


v  ■  2xc[f(r)cos$  *  g(r)s1n$] 


f"  ♦  *7  f ^7  f  -  Re  g  ■  -  Re  r 

9"+79'  -  779  +  Ref*0 

#_.-0  at  r  »  1 

T’g  1  finite  at  r  -  0. 


(4) 

(5a) 

(5b) 

(5c) 


The  prime  denotes  d/dr.  These  equations  seem  to  call  for  an  analytical 
solution,  but  our  attempts  were  not  yet  successful.  For  low  Reynolds  numbers. 
It  is  obvious  that 
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f  -  n(Rnl,  n  -  n(Ro?)  (6) 

and  In  more  detail 

f  •  ^2.  (r  -  r3)  +  0(Re3)  (7a) 

9  "  if?  <2r  *  3r3  *  pS)  +  °(Re4)-  (7b) 

A  series  Including  higher  order  terms  In  Re  has  been  constructed  but  converges 
only  for  Re<12.  Another  Interesting  limit  can  be  easily  obtained  for  Re-*—: 


f9lr  )  (8) 

Owing  to  the  loss  of  the  highest  derivatives,  this  solution  cannot  satisfy  the 
boundary  conditions  and  Is  valid  only  outside  the  thin  boundary  layers  near 
the  cylinder  sidewall. 

Without  any  detailed  knowledge  of  the  solution  In  the  medium  range  of  Re, 
It  can  be  seen  that  the  basic  flow  exhibits  characteristically  different 
behavior  at  low  and  high  Reynolds  numbers.  At  low  Re,  the  component  f  in  the 
plane  <,»0  (Z.z-plane)  Is  dominating.  At  high  Re,  f  is  negligible  (except  In 
the  boundary  layer),  and  the  component  g  In  the  plane  $-90°  Is  dominating. 

In  the  medium  range  of  Re,  we  have  applied  a  spectral  collocation  method 
for  numerical  solution.  This  method  is  especially  capable  of  resolving  the 
steep  gradients  in  the  boundary  layers  at  large  Re.  Figure  1  shows  the  numer¬ 
ical  results  In  the  form  of  contour  lines  of  equal  axial  velocity.  The  flow 
to  the  right  or  above  the  (thicker)  zero-velocity  line  Is  toward  the  reader. 
The  shift  of  the  velocity  maximum  from  $*0  at  Re»0.01  to  $«90°  at  R-1000,  and 
the  pronounced  boundary  layer  formation  In  the  latter  case  Is  clearly  visible. 


Figure  1.  Contour  lines  of  equal  axial  velocity,  vz/(2to)  «  const.,  for 
Re-0.01,  10  and  1000.  The  difference  between  levels  Is  0.001,  0.1  and  0.2, 
respectively. 
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for  supporting  the  relevance  of  our  solution.  Figure  2  shows  the  compari¬ 
son  with  computational  results  for  the  velocity  distrihut  ions  in  t.ho  center 
cross-section  (2=0)  of  the  cylinder.  The  data  were  kindly  provided  by  Dr. 
Vaughn,  Sandia  National  Laboratories.  The  agreement  for  Re=14.9  is  considered 
representative  for  the  range  of  lower  Reynolds  numbers.  At  the  higher  Re=45.7, 
a  minor  but  systematic  deviation  between  our  solution  and  the  Sandia  results 
exists  that  seems  to  be  due  to  a  superposed  cellular  motion.  This  cellular 
motion  is  not  yet  incorporated  Into  In  our  analysis  and  will  be  subject  to 
further  study.  The  axial  flow  considered  here  will  provide  the  basic  flow  for 
an  analysis  of  the  hydrodynamic  stability.  It  is  encouraging,  though,  that 
our  simple  theory  yields  results  in  essential  agreement  with  the  computational 
solution  of  the  full  Navler-Stokes  equations  for  a  finite-length  cylinder.  The 
computational  effort  for  solving  the  ordinary  differential  equations  (5)  is 
rather  small,  typically  25  msec  (IBM  3081)  per  solution. 


Figure  2.  Radial  distribution  of  the  axial  velocity  (in  fps)  at  2=0  for 
Re»14.9  and  Re*>45.7.  The  symbols  show  the  solution  to  the  full  Navler-Stokes 
equations  (Vaughn,  personal  communication). 


Moments 

In  order  to  compare  with  the  experimental  data  for  the  despin  moments,  we 
have  analyzed  the  moments  that  arise  from  the  deviation  velocity  Vd,  disre¬ 
garding  the  contribution  of  pure  solid  body  rotation.  Conservation  of  angular 
momentum  for  the  fixed  cylindrical  control  volume  R  (surface  S)  requires 


5  -  It  k  •  Li)odR  *  JJJ  k 


(2 2  *  Vd) ]pdR 


+ 


(L  +  v.)]p(v 


n)dS 


(9) 
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where  M  Is  the  resultant  torque  on  the  system.  The  first  term  on  the  right- 
hand  side  vanishes  for  steady  The  second  term  originates  from  Coriolis 
forces  In  the  nutating  system.  The  third  term  describes  the  net  rate  of 
angular  momentum  flux  through  the  control  surface  with  the  outer  normal  n. 
The  contribution  from  £  *  ^  vanishes  since  ^  has  only  an  axial  component. 
The  last  term  gives  the  moments  due  to  the  shear  forces  at  the  cylindrical 
surface.  There  Is  no  pressure  contribution  since  pd»0. 

Substitution  of  and  v2  provides  the  following  expressions  for  the 
components  of  M: 


mJt(2Q«s1neHuB)  -^J1  r2fdr} 

(10a) 

m1(2has1ne)(«a)  (  -  ~^1^-  -pJ  1  r2gdr} 

(10b) 

m,(2nas1ne)2  {  J 1  r2fdr} 

1  0 

(10c) 

where  m,  Is  the  liquid  mass  In  the  cylinder.  The  derivatives  g',  f*  In  M^  My 
originate  from  the  shear  forces  at  the  cylinder,  whereas  the  Integra's  provide 
the  rate  of  angular  momentum  flux.  The  despin  moment  M  Is  due  to  Coriolis 
forces.  Note,  however,  that  the  Interpretation  of  Mz  will  be  different  If  the 
flux  through  the  cylinder  end  walls  Is  prohibited. 


Using  the  differential  equations  (5a, b)  and  integration  by  parts,  the 
Integrals  In  (lOa.b.c)  can  be  expressed  as 


r2fdr  -  -  Sq.U) 


•,  r2gdr 


T'(l)  +  * 
“Re  +  ** 


(ID 


The  moments  can  therefore  be  written  as 


mt(2cas1ne)(wa)  M  , 

(12a) 

mi(2nas1ne)2  Mz,  Mz 

•-is'd) 

(12b) 

(i)  ^ 

nslne  z* 

(120 

Since  f * ( 1 1  <  0,  g* ( 1 )  <  0  for  all  finite  Re<0,  the  moment  Mz  about  the  spin 
axis  Is  always  positive,  whereas  Mx  <  0.  For  small  Reynolas  number,  the 
series  solution  (7)  provides 
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which  allows. a  quick  estimate  of  the  moments  for  Re  <  10,  say.  The  linear 
Increase  of  M2  with  Re  Is  consistent  with  the  experimental  data. 

The  values  of  My  and  hL  frg«  numerical  solutions  are  given  In  Figyre  3. 
It  Is  Interesting  tcrnote  tnat  M2  has  a  pronounced  maximum  at  Re-19.  Mv 
varies  between  1/4  and  -1/4  with  a  sign  change  at  Re-31.  A  comparison  tilth 
experimental  data  for  K_  (Miller  1982)  Is  glyen  in  Figure  4  on  a  doubly 
logarithmic  scale.  Thezdata  are  reduced  to  M,  using  the  initial  spin  rate 
(0-4000  rpffl.  For  Re  <  10,  the  experimental  points  agree  with  the  simple  law 
M*  ■  Re/96.  The  systematic  deviation  for  higher  Reynolds  numbers  Re  <  200 
seems  to  be  due  to  the  occurrence  of  a  cellular  motion  that  Is  likely  to 
reduce  the  moment  M-.  The  two  data  points  for  Re  >  103  Indicate  turbulent 
flow.  2 


Figure  3.  The  nondlmenslonal  coefficients  tL  and  M,  In  the  moments  (12a.b) 
versus  the  Reynolds  number  Re.  J 


Figure  4.  Comparison  gf  the 
theoretical  values  of  Mz 
with  the  experimental  data  of 
Miller  (1982).  The  straight 
}1ne  shows  the  asymptotic  law 
M7  -  Re/96. 


Conclusions 

The  theory  developed  here  Is  a  first  but  virtually  essential  step  toward 
understanding  and  predicting  the  gross  features  of  the  fluid  motion  in  a 
spinning  and  nutating  cylinder.  The  axial  basic  flow  allows  an  analysis 
of  the  hydrodynamic  Instability  with  respect  to  a  cellular  motion  and  provides 
the  basis  for  refinements  of  the  theory. 
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HIGHLY  VISCOUS  FLUID  FLOW 
IN  A  SPINNING  AND  NUTATING  CYLINDER 


Thorwald  Herbert 

Department  of  Engineering  Science  and  Mechanics 
Virginia  Polytechnic  Institute  and  State  University 
81ack$burg,  Virginia  24061 

ABSTRACT.  Spln-stabll 1 2ed  projectiles  with  liquid  payloads  can  experi¬ 
ence  a  severe  flight  instability  characterized  by  a  rapid  yaw  angle  growth  and 
a  simultaneous  loss  In  spin  rate.  Laboratory  experiments  and  field  tests  have 
shown  that  this  instability  originates  from  the  Internal  fluid  motion  In  the 
range  of  high  viscosity.  Evaluation  of  the  experimental  data  and  analysis  of 
the  equations  for  the  fluid  motion  In  a  spinning  and  nutating  cylinder  suggest 
a  theoretical  approach  In  three  major  steps:  (1)  analysis  of  the  steady  vis¬ 
cous  flow  In  an  Infinitely  long  cylinder,  (2)  hydrodynamic  stability  analysis 
of  this  basic  flow,  and  (3)  analysis  of  the  end  effects.  The  basic  flow  has 
been  found  in  analytical  form.  At  low  Reynolds  number,  this  flow  agrees  well 
with  computational  results  for  the  center  section  of  a  cylinder  of  aspect 
ratio  4.3.  The  despin  moment  caused  by  this  flow  largely  agrees  with  experi¬ 
mental  data  for  a  wide  range  of  Reynolds  numbers.  Current  work  alms  at  the 
stability  of  this  flow. 

1.  INTRODUCTION .  It  is  well-known  that  spin-stabilized  shells  carrying 
liquid  payloads  can  suffer  dynamical  instability.  For  cylindrical  cavities 
and  low  viscosity  of  the  liquid,  the  Instability  due  to  basically  Invlscld 
Inertial  waves  can  be  predicted  by  the  Stewart son-Wedemeyer  theory  (1,21. 
This  theory  rests  on  the  boundary- layer  approach  and  is,  therefore,  restricted 
to  the  range  of  sufficiently  large  Reynolds  numbers.  The  Instability  of  cer¬ 
tain  shells  like  the  XM  761  [3,4],  howeve-*,  escaped  such  a  prediction  and  Is 
also  distinguished  In  character  owing  to  the  rapid  loss  In  spin  rate.  Experi¬ 
ments  with  a  full-scale  liquid  cylinder  (5)  and  subsequent  field  tests  (6] 
established  that  this  new  flight  Instability  Is  most  pronounced  for  liquid- 
fills  of  very  high  viscosity. 

We  conduct  a  theoretical  analysis  of  this  problem  In  order  to  support  the 
ongoing  experiments  and  to  independently  obtain  Insight  Into  the  anatomy  of 
the  flow  phenomena.  The  Inttial  steps  of  this  analysis  are  reported  elsewhere 
(7):  evaluation  of  the  experimental  data  base,  dimensional  analysis,  scaling 
aspects,  governing  equations,  and  discussion  of  various  simplifying  assump¬ 
tions.  Two  observations  in  this  earlier  work  led  to  the  building-block 
approach  discussed  In  the  following.  First,  if  the  despin  (negative  roll) 
moments  (5)  and  void  observations  18]  are  correlated  with  the  Reynolds  number 
Re,  at  least  three  regions  can  be  distinguished.  At  low  Re,  the  despin  moment 
increases  proportional  to  Re,  and  the  void  In  an  incompletely  filled  cylinder 
Is  parallel  to  the  spin  axis.  This  suggests  a  simple  fluid  motion  that  Is 
essentially  independent  of  the  axial  coordinate,  except  In  the  neighborhood  of 
the  end  walls.  In  a  middle  range  of  Re,  the  despin  moment  assumes  a  maximum, 
and  a  wavy  distortion  of  the  void  seems  to  indicate  a  cellular  structure  of 
the  fluid  motion.  This  cellular  motion  can,  In  principle,  originate  from 
hydrodynamic  instability  of  the  basic  flow  with  respect  to  axially  periodic 
disturbances.  At  still  higher  Reynolds  numbers,  the  despin  moment  decreases 
with  increasing  Re  In  a  manner  not  clearly  defined  by  the  few  available  data 
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points.  The  void  observations  indicate,  however,  that  the  motion  ultimately 
becomes  turbulent. 

The  second  observation  was  the  appearance  of  the  nutation  rate  and  angle 
as  a  small  parameter  In  the  equations  for  the  deviation  from  solid-body  rota¬ 
tion.  The  forcing  term  due  to  nutation  can  be  considered  small  enough  for 
linearization  of  the  equations  In  the  situations  of  practical  and  theoretical 
Interest. 

In  the  following,  we  describe  the  development  of  a  simple  system  of  equa¬ 
tions  for  the  basic  flow.  Analytical  solutions  are  given  for  the  flow  field 
and  for  the  liquid  moments.  A  comparison  Is  made  with  computer  simulations  of 
the  flow  (91  and  with  experimental  data  (51.  The  properties  of  Inertial  modes 
at  low  Reynolds  numbers  and  the  possibility  of  Instability  due  to  primary 
resonance  is  discussed. 

2.  GOVERNING  EQUATIONS.  We 
consider  the  motion  of  a  fluid  of 
density  0  and  viscosity  y  In  a  cyl¬ 
inder  of  radius  a  and  length  2c  that 
rotates  with  the  spin  rate  u  about 
its  axis  of  symmetry,  the  z-axls. 

We  consider  the  motion  with  respect 
to  the  nutating  coordinate  system 
x.y.z.  This  system  Is  obtained  from 
the  Inertial  system  X.Y.Z  by  a  rota¬ 
tion  with  the  nutation  angle  e  about 
the  axis  Y-y.  Therefore,  x  Is  In 
tne  Z,z-planef  and  this  plane 
rotates  about  the  Z-axls  with  the 
nutation  rate  a.  The  two  axes  of 
rotation  Intersect  In  the  center  of 
mass  of  the  cylinder,  as  shown  in 
Fig.  1.  In  contrast  to  the  experi¬ 
mental  procedures  (51,  we  consider 
w>0,  q,  and  0s8<*/2  as  constant.  The 
fluid  motion  is  governed  by  the 
Navler-Stokes  equations  written  In 
the  nutating  coordinate  system: 

°V 

°lot  ■  20  *  yn  +  2  *  (2  *  Ol  •  -vPn  +  uv*vn 

v.yn  -  o. 

Vn  Is  the  velocity  measured  in  the  nutating  frame,  P  the  pressure,  and  r  the 
position  vector.  The  body  force  due  to  gravity  has  been  disregarded.  Equa¬ 
tions  (1)  are  subject  to  the  no-slip  and  no-penetration  conditions  at  the 
cylinder  walls. 

It  Is  convenient  (7)  to  split  the  velocity  and  pressure  fields  according 
to 


4-  Z 


spin  nutation 
axis  axis 


Figure  1.  Definition  sketch. 
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where  ys,P_  describe  the  state  of  pure  solid  body  rotation,  whereas 
represent  the  deviation  from  solid  body  rotation.  The  advantage  of  this  iso¬ 
lated  view  on  the  deviation  is  obvious:  Vj  and  the  reduced  pressure  Pd  are 
responsible  for  the  observed  flight  instability.  A  glance  at  the  equations 
shows  that  Vd=0  and  Pd=0  If  either  one  of  the  following  conditions  Is  satis¬ 
fied:  w*=0,  u«0,  e-0  or  u*®  (solid  fill). 

The  equations  for  Vd,  Pd  are  then  written  In  terms  of  nondlmenslonal 
quantities  vd,  pd.  We  use  a,  u,  and  o  for  scaling  length,  time,  and  mass. 
Note  that  this  choice  Is  ambiguous  17)  and  excludes  the  case  u»0  (which  lacks 
practical  interest).  The  problem  then  depends  on  four  nondlmenslonal 
parameters: 


x  =  c/a 
a  *  sine 
t  *  n/w 

Re  *  pwaJ/u 


aspect  ratio 

frequency 
Reynolds  number. 


The  aspect  ratio  e..ters  the  solution  only  through  the  boundary  conditions.  The 
boundary  conditions  on  yd  are  homogeneous. 

In  cylindrical  coordinates  r,o,2,  the  equations  for  the  nondlmenslonal 
deviation  velocity  yd  =  (vr,v  ,v2)  and  pressure  pd  take  the  form 

i  ,  .  !v.  3v 

7F<rvr>  +f  ^  +  =  0  <2a> 

D'v  -  -r  -  2(1  +  x  )  v  +  2t  v 

r  r  za  oz  (2b) 

3prt  I  Vr  ?  3V* 

* - 1  +  i-  [0"v  -  —  -  - - £) 

ar  Re  r  r2  r2  a$ 


where 


O'v  + 


+  2(1  +  ^z)  vr  -  2irvz 


1  3PH  1  V*  ?  3Vr 

a  .  I  _!  +  1  I  0"V  .  -i  +  i - H] 

¥•  Si  Da  1  a  si  * 


D'v  +  2t  v  -  2t  v 


'  ri  r 2 
3pd  1 

JT  -  2rv +  k  D"va 


D'  =  i—  +  i—  +  v  L  +  .1  L 

at  a$  r  ar  r  a$ 


o-  =  +  1 1-  ♦  li¬ 

ar*  r  3r  r*  a** 
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and 


t  =  -tcos$,  t  *  es1n$,  t,  =  icose,  c  =  tsine  (3) 
r  *  z 

The  primary  effect  of  nutation  is  contained  in  the  ♦-periodic  force  term 
-2rir  *  2crcos$  In  the  z-momentum  equation  (2d).  If  this  term  vanishes 
throughout,  t-0,  equations  (2)  support  a  trivial  solution  yd“0,  pd*0.  For 
sufficiently  small  e*Q,  It  Is  obvious  that  the  deviation  velocity  Is  of  order 
0(c).  In  the  situations  of  practical  interest,  c  --  (o/w)s1ne  turns  out  to  be 
a  rather  small  parameter.  Even  a  conservative  estimate  with  o  s  500  rpm, 
u  i  3000  rpm,  and  a  s  20°  provides  values  of  c  s  0.057.  Consequently,  it 
seems  well  justified  to  linearize  the  equations  In  e.  This  linearization 
Imposes  no  restriction  on  the  Reynolds  number. 

3.  THE  BASIC  FLOW.  The  system  of  equations  after  linearization  Is  still 
quite  difficult  to  solve.  Any  serious  attempt  to  satisfy  all  boundary  condi¬ 
tions  leads  directly  to  a  purely  computational  approach.  Use  of  the  boundary- 
layer  approximation  would  simplify  the  task  but  seems  Inappropriate  In  the 
Interesting  range  of  low  Reynolds  numbers.  Recalling  that  the  flow  In  a  rela¬ 
tively  long  cylinder  (aspect  ratio  x«4.3)  at  low  Re  exhibits  little  axial 
variation  over  much  of  the  cylinder  length  [7],  we  have  relaxed  the  boundary 
conditions  at  the  end  walls.  As  a  first  step,  we  seek  for  a  steady  flow  In  an 
Infinitely  long  cylinder. 

At  closer  analysis,  the  z-lndependent  force  term  In  eq.  (2d)  can  only  be 
balanced  by  a  purely  axial  deviation  velocity.  It  Is  consistent  with  the  lin¬ 
earized  equations  to  assume  a  solution  In  the  form 


yd  -  (o,  o,  vz),  Pd  -  o 

(4) 

and  moreover. 

v2  "  vz<r’*>  "  Mf(r)cos*  +  g(r)slna) 

(5) 

Substituting  (4), (5)  into  the  linearized  equations  provides 

f"  +  7  f’-  77  f  -  Re  g  •  -  Re  r 

(6a) 

g"  tig'  -  f,  g  ♦  Re  f  •  0 

(6b) 

f=0,  j-0  at  r=l 

(6c) 

f,  g  finite  at  r=0 

(6d) 

The  primes  denote  d/dr.  For  Re*0,  the  solution  of  these  equations  can  be 
found  In  the  form  of  series 


f  -  5s  (r  -  r3)  -  (7r  -  12r°  +  6r*  -  r’)  +  0(Re*)  (7a) 
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9  =  fff  {2f  -  3r J  +  rs)  +  O(Re-).  (7b) 

With  higher  terms  included,  these  series  converge  for  Re  <  12.  In  the  limit 
Re  -  «>,  one  obtains 


f  •  0,  g  «  r  •  as  Re  ♦  (8) 


Owing  to  the  loss  of  the  highest  derivatives,  however,  this  solution  cannot 
satisfy  the  boundary  conditions  (6c)  and  is  valid  only  outside  the  thin  bound¬ 
ary  layers  near  the  wall  at  r=l.  Even  without  any  knowledge  of  the  solution 
in  the  intermediate  range,  the  different  character  of  the  basic  flow  at  low 
and  high  Reynolds  numbers  Is  evident.  At  low  Re,  the  component  f  in  the  z,x- 
plane  <t>a0  Is  dominating.  At  high  Re,  f  Is  negligible  In  the  core  of  the 
cylinder  while  g  In  the  z,y-p1ane  $=90°  is  dominating. 

In  earlier  work  (10),  we  have  applied  a  spectral  collocation  method  for 
numerical  solution  of  eqs.  (6).  Here,  we  derive  an  analytical  solution  by 
introducing  the  complex  function  F=g+1f.  Eqs.  (6)  can  then  be  written  In  the 
form 


+  rF'  -  (1  +  1Re  r*)F  =  -  iRe  r3 

(9a) 

F  =  0  at  r  =  1 

(9b) 

F  finite  at  r  *>  0 

(9c) 

A  particular  solution  of  the  Inhomogeneous  equation  (9a)  Is  F0=r,  whereas  the 
homogeneous  part  of  (9a)  Is  the  equation  for  the  modified  Bessel  functions 
1 1 (qr )  and  Kt(qr)  of  the  complex  argument  qr  with  q  *  /Re/2  (1  +  1)r.  For 
(9c),  Kj  (qr)  cannot  contribute  to  the  solution.  Finally,  (9b)  provides 

F(r)  *  g  4  if  *  r  -  I t (qr)/I ^ (q) .  (10) 

Expressing  the  solution  in  terms  of  Kelvin  functions  of  real  argument  Is  of 
little  advantage  for  the  numerical  evaluation.  The  solution  Is  valid  for 
arbitrary  Re  but  may  be  unstable  as  Re  exceeds  some  critical  value.  It  Is 
straightforward  to  derive  the  approximations  (7)  from  the  ascending  series  for 
I.  (and  to  explain  the  convergence  problem  for  larger  Re).  The  asymptotic 
expansion  for  large  arguments  provides 

F  -  r  -  /r  eq(r_1).  (11) 

This  expression  agrees  to  within  with  (8)  provided  that  r  <  1  -  4.  The 
boundary  layer  thickness  4  can  be  obtained  from  the  transcendental  equation 

6  =  /2/Re  [4,605  -  |  an  (1  -  4)),  (12) 

e.g.,  4=0.223  for  Re=1000.  The  characteristic  changes  in  the  flow  structure 
with  increasing  Re,  In  particular  the  shift  of  the  velocity  maximum  from  $»0 
at  Re=2  to  $*90°  at  Re=200  are  shown  In  Fig.  2. 
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Fig.  3  compares  the  dimensional  velocity  distributions  obtained  from  (5), 
(10)  with  computational  results  for  the  center  cross-section  (z=0)  of  the 
cylinder.*  The  agreement  for  Re=14.9  Is  considered  representative  for  the 
range  of  lower  Reynolds  numbers.  The  numerical  simulation  provides  very  small 
components  vr,  v  and  hence  verifies  our  estimates.  At  the  higher  value 
Re*45.7,  a  systematic  deviation  between  the  two  results  seems  to  be  due  to  a 
superposed  cellular  motion  that  is  not  yet  Incorporated  Into  our  analysis.  It 
Is  encouraging,  however,  that  the  simple  theory  of  the  basic  flow  yields 
results  In  essential  agreement  with  the  computational  solution  of  the  full 
Navler  Stokes  equations  for  a  finite-length  cylinder. 


4.  MOMENTS.  With  the  deviation  velocity  Vj  -  (0,0,wav2)  and  vz  given, 
the  moments  on  a  finite-length  section  of  the  cylinder  can  be  calculated.  We 
consider  a  control  volume  R  (surface  S)  formed  by  the  solid  cylindrical  wall 
and  liquid  surfaces  at  both  ends.  Conservation  of  angular  momentum  requires 

M  +  t  (r  «  Fd)  -  ft  JJ7 (r  *  yd)odR  +  JJJ  (r  «  (2a  ,  yd)]0dR 

R  R 

(13) 

+  J.r  (r  *  yd)o(Vd  •  n)dS  ♦  JJ  (r  „  ys)o(Vd  •  n)dS 

where  n  Is  the  outer  unit  normal.  On  the  left-hand  side,  M  Is  the  resultant 
torque'on  the  control  volume.  The  second  term  accounts  for  the  moments  due  to 
the  shear  force  Fd  and  vanishes  owing  to  the  solid  sidewall  and  cancellation 
of  the  contributions  from  both  ends**.  On  the  right-hand  side,  the  first  term 
vanishes  for  steady  Vd.  The  second  term  originates  from  Coriolis  forces  in  the 
nutating  system.  The  third  term  vanishes  since  yd  has  only  an  axial  compo¬ 
nent.  The  last  term  then  provides  the  net  rate  of  angular  momentum  flux 
through  the  control  surface. 

Substitution  of  Vd  leads  to  the  following  expressions  for  the  cartesian 
components  of  M: 

M  ■  m  (2na$1ne)(wa)  m  ,  m  »  -  /'  rJfdr  (14a) 

A  1C  n  A  q 

Mu  =  m  (2nas1ne)(uja)  m  ,  m  -  ■  r*gdr  (14b) 

y  i  m  y  y  o 

M  »  m  (2nas1ne)'  m  ,  m  =  /'  r^fdr  «  -m  (14c) 

■  <C  Z  Z  Q 

where  m^  Is  the  liquid  mass  In  the  cylinder.  In  this  form,  the  components 
Mx,  M  represent  the  net  rate  of  angular  momentum  flux  through  the  liquid  end- 
walls,  whereas  the  roll  moment  Mz  Is  solely  due  to  Coriolis  forces.  A  close 


*The  data  were  kindly  provided  by  Or.  H.  Vaughn,  Sandla  National  Laboratories. 

♦♦Improper  account  of  the  sidewall  conditions  Introduced  an  Incorrect  factor 
of  two  In  earlier  results  for  Mx,  ( 10 | . 
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relation  between  roll  moment  My  and  yaw  moment  Mx  has  also  been  found  by 
Murphy  (11)  for  the  range  of  hign  Reynolds  numbers. 

A  different  Interpretation  can  be  derived  using  the  differential  equation 
(9a),  Integrating  by  parts,  applying  (9b),  and  separating  real  and  i..iag1nary 
part: 

mz  "  ‘mx  =  ^2fdr  *  "  (15a) 

«y  "  -  Jl  r’gdr  -  -  (15b) 


In  this  form,  the  moments  are  directly  related  to  the  shear  forces  at  the 
cylindrical  sidewall,  r«l.  Since  f'(l)<0,  g'(l)<0,  the  roll  moment  Mz  Is 
always  positive  (even  for  n<0),  while  Mx  Is  negative  for  n<0  and  changes  Sign 
with  o.  For  small  Re,  the  series  (7)  provide  the  approximations 


_  Rt.  _  Rej 

z  “  95’  y  "  “  T535 


that  can  be  used  for  quick  estimates  up  to  ReslO.  The  linear  Increase  of  mz 
and  M  with  Re  Is  consistent  with  the  experimental  data.  From  the  analytical 
solution  (10),  we  obtain 


F'(l)  -  g ' ( 1 )  +  if'  (1)  -  2  -  ql0(q)/l1(q).  (17) 

Substitution  into  (15)  provides  the  variation  of  mz,  m  with  the  Reynolds 
number  shown  In  Fig.  4.  The  coefficient  mz  assumes  a  pronounced  maximum  at 


Fig.  4.  The  nondimenslonal  coef-  Fig.  5.  Comparison  of  the  theoretl- 
fldents  m  ,  m  In  eq.  (14)  versus  cal  result  for  m  with  experimental 
the  Reynolds  number,  Re.  data  (5].  The  straight  line  snows 

the  asymptotic  law  m  »Re/96. 
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Re«19.  The  coefficient  itiy  Is  negative  and  reaches  an  asymptotic  value  of 
my  -  -  1/4  as  Re  -  «.  Renee,  for  n>0,  My  reduces  the  pitch  moment  due  to 
the  solid  body  rotation.  This  result  Is  consistent  with  computations  19).  A 
comparison  of  theoretical  and  experimental  (5)  results  for  m.  Is  given  In  Fig. 
5  on  a  double  logarithmic  scale.  The  Initial  spin  rate  u*4000  rpm  was  used 
for  reducing  the  experimental  data.  For  Re<10,  the  experimental  points  match 
the  analytical  result  as  well  as  the  asymptotic  law  m,  *  Re/96.  The 
occurrence  of  a  maximum  of  mz  Is  found  to  be  a  property  of  the  basic  flow. 
Only  the  systematic  deviation  for  higher  Reynolds  numbers  Re<200  may  be 
attributed  to  a  cellular  motion.  The  two  data  points  for  Re>10*  probably 
Indicate  turbulent  flow. 

The  basic  flow,  hence,  can  be  considered  a  first  but  essential  step 
toward  understanding  and  predicting  the  gross  features  of  the  fluid  motion  In 
a  spinning  and  nutating  cylinder.  Some  observations,  however,  such  as  the 
virtual  Independence  of  the  despin  moment  on  the  spin  rate  require  further 
analysis,  especially  of  the  end  effects.  The  occurrence  of  a  cellular  motion 
may  be  due  to  hydrodynamic  Instability  of  the  basic  flow. 

5.  STABILITY  ANALYSIS.  The  stability  analysis  Is  currently  coi-ucted 
and  only  a  brief  outline  Is  given  here.  We  superpose  to  the  steady  flow 
vn*(o,r,vz),  pn  disturbances  v'*(u,v,w),  p  sufficiently  small  for  lineariza¬ 
tion.  Substitution  Into  eqs."(2)  and  neglect  of  products  between  disturbances 
and  terms  of  order  0(eJ)  provides  the  following  stability  equations: 


‘It  *  w  *  lf>  -  <°-u  -  77  -  77 

*  (,2  «  *  V  ■  0 

*  I? +  2<‘  *  'z)u  *  ?  H>  -  h  ^  *  77  -  fr  tl> 

‘"'ll-  2'r“>  ■  0 

*  I!  ♦  •  h  <0■“,  +  -  2VU  +  <7  TT  *  2,r>v 

*  vz  S>  *  0 


(18®) 


(18b) 


(18c) 


I  L.  (ru)  ♦  I  2*  ♦  w  .  o 

r  3r  '  '  r  36  3z 


(18d) 


Three  groups  of  terms  have  been  separated  by  braces  in  eqs.  (18a-c).  The 
first  group,  if  set  to  zero,  represents  the  equations  for  Invlscid  Inertial 
modes  ~exp(1m<t>  +  laz  +  st),  where  m  Is  the  (Integer)  azimuthal,  a  the  axial 
wavenumber,  and  s  =  sr  +  1s^  provides  the  amplification  rate  sr  (»0)  and  fre¬ 
quency  $i.  Usually,  an  equation  for  the  pressure  Is  used  for  obtaining  the 
analytical  solution.  We  have  derived  an  alternative  system  in  terms  of  u,v 
and  applied  the  spectral  method  to  be  used  for  more  general  cases  In  order  to 
check  the  numerical  results  against  the  exact  values. 
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The  second  group  of  terms  multiplied  with  1/Re  represents  the  viscous 
correction  to  the  Inertial  modes.  By  eliminating  w  and  p,  a  system  of 
ordinary  differential  equations  for  u  and  v  has  been  derived.  Due  to  the 
higher  order  of  this  system,  all  boundary  conditions  can  be  satisfied.  Pro¬ 
grams  have  been  developed  for  calculating  spectra  of  (complex)  eigenvalues, 
for  tracing  single  eigenvalues  as  function  of  Re,  a,  and  for  obtaining  the 
eigenfunctions.  At  high  Re,  the  results  follow  the  trends  predicted  by 
asymptotic  theories.  Our  analysis,  however,  also  covers  the  range  of  low 
Reynolds  numbers,  where  the  Inertial  modes  suffer  rapid  decay  (sr<0). 

The  most  Interesting  aspect  of  the  stability  equations  (18)  Is  the  third 
group  of  terms.  The  coefficients  In  this  group,  vz,  tr,  and  t*»  are  (1)  of 
order  0(c)  and  (2)  periodic  In  a.  The  periodicity  In  a  leads  to  a  coupling  of 
the  mode  equations  for  m  and  mil,  and  may  cause  primary  resonance  between 
Inertial  modes.  In  view  of  viscous  damping,  this  resonance  Is  likely  to  occur 
as  c  exceeds  a  critical  value  that  decreases  as  Re  Increases.  The  analysis  of 
this  parametric  Instability  Is  currently  In  preparation. 
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Abstract 

Spin-stabilised  projectiles  with  liquid  payloads  can  experience  a  severe  flight  instability  character¬ 
ised  by  a  rapid  yaw  angle  growth  and  a  simultaneous  loss  in  spin  rate.  Laboratory  experiments  and  Geld 
tests  have  shown  that  this  instability  originates  from  the  internal  fluid  motion  in  the  range  of  high 
viscosity.  We  have  developed  a  simple  model  of  this  internal  motion  that  provides  the  flow  field  and  the 
liquid  moments  in  ansdytical  form.  A  detailed  comparison  of  the  roll  moment  with  data  from  spin-down 
equipments  is  given.  New  experimental  procedures  are  suggested. 

1.  Introduction 

It  is  well-known  that  spin-atabilited  shells  carrying  liquid  payloads  can  suffer  dynamical  instability. 
For  cylindrical  cavities  and  low  viscosity  of  the  liquid,  the  instability  due  to  basically  inviscid  inertial 
waves  can  be  predicted  by  the  Stewartson-Wedemeyer  theory.  The  instability  of  certain  shells  like  the 
XM761,  however,  escapes  such  a  prediction  and  is  also  distinguished  in  character  by  the  rapid  loss  in 
spin  rate.  Experiments  with  a  full-scale  liquid-filled  cylinder  [l]  and  subsequent  field  tests  [2]  establish 
that  this  new  flight  instability  is  most  pronounced  for  liquid  fills  of  very  high  viscosity. 

We  have  conducted  s  theoretical  analysis  of  the  problem  with  special  attention  to  the  range  of 
high  viscosity.  Two  observations  have  permitted  the  development  of  a  simple  model  of  the  internal 
fluid  motion.  First,  the  nutation  rate  and  angle  appear  as  a  small  parameter  in  the  nondimensional 
equations  for  the  deviation  from  solid-body  rotation  ( 3] .  Therefore,  the  forcing  term  due  to  nutation 
can  be  considered  sufficiently  small  for  linearisation  in  the  situation  of  practical  interest.  Second,  in  a 
sufficiently  long  cylinder,  the  velocity  field  consists  of  primarily  an  axial  component  that  is  independent 
of  the  axial  position  (4).  The  other  components  are  of  the  same  order  only  in  the  neighborhood  of  the 
end  walls.  Therefore,  essentia]  features  of  the  internal  fluid  motion  can  be  obtained  by  studying  the 
flow  in  a  finite  segment  of  an  infinitely  long  cylinder. 

Although  the  agreement  of  the  theoretical  and  experimental  results  on  the  roll  moment  is  surpris¬ 
ingly  good,  there  seems  to  be  an  essential  discrepancy:  while  the  theoretical  result  depends  (via  the 
Reynolds  number)  on  the  spin  rate,  Miller  (lj  found  that  "the  despin  moment  was  not  a  function  of  the 
canister  spin  rate,  provided  a  sufficient  spin  rate  ia  present.  With  more  detailed  data  available  in  a  wide 
range  of  Reynolds  numbers,  we  shed  some  light  on  this  virtual  disagreement.  The  experiments'  data 
follow  in  fact  the  theoretically  predicted  trends,  A  systematic  deviation  persists  in  the  range  of 
extremely  high  viscosities,  This  deviation  is  likely  to  originate  from  changes  in  temperatures  and  con¬ 
sequently  in  viscosity.  It  also  turns  out  that  a  survey  of  the  roll-moment  versus  Reynolds  number 
curve  can  be  obtained  in  very  few  experimental  runs  at  different  viscosity. 
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2.  Analysis  at  the  Flow  Field 

We  consider  the  motion  of  s  fluid  of  density  p  and  viscosity  p  in  a  cylinder  of  radius  a  and  length 
2e  that  rotates  with  the  spin  rate  w  about  its  axis  of  symmetry,  the  r-axis.  We  consider  the  motion 
with  respect  to  the  nutating  coordinate  system  x,y ,z.  This  system  is  obtained  from  the  inertial  system 
X,Y,Z,  by  a  rotation  with  the  nutation  angle  8  about  the  axis  Y=y.  Therefore,  x  is  in  the  Z,*-plane, 
and  this  plane  rotates  about  the  Z-axis  with  the  nutation  rate  fl .  The  two  axes  of  rotation  intersect  in 
the  center  of  mass  of  the  cylinder,  as  shown  in  Fig.  1.  In  contrast  to  the  experimental  procedures  [l], 
we  consider  w  >0,  fl  ,  and  0 <8  ^tr  /2  as  constant. 

It  is  convenient  [3]  to  split  the  velocity  and  pressure  fields  according  to 

V.  =  V.  +  V,t  P.=  P,  +  Pi  (1) 

where  V,  ,P,  describe  the  state  of  pure  solid  body  rotation,  whereas  \i,Pi  represent  the  deviation  from 
solid  body  rotation.  The  advantage  of  this  isolated  view  on  the  deviation  is  obvious:  and  the  reduced 

pressure  Pi  are  responsible  for  the  observed  flight  instability.  A  glance  at  the  equations  [3]  shows  that 
Vi  ss  0  and  Pi  k  0  if  either  ui—0,  fl  =0,  0  =  9,  or  p  — *oo  (solid  fill). 

The  effect  of  nutation  and  hence  the  deviation  velocity  is  of  order  0(e)  where  <  =  (fl  /w)sin0.  In 
the  situations  of  practical  interest,  <  turns  out  to  be  rather  small:  a  conservative  estimate  provides 
values  of  <<0.054.  Consequently,  it  seems  well  justified  to  linearise  the  equations  in  e.  This  linearisa¬ 
tion  imposes  no  restriction  on  the  Reynolds  number  Re  =  pua2/p  .  Recalling  that  the  flow  in  a  rela¬ 
tively  long  cylinder  (aspect  ratio  X  =  4.3)  at  low  Re  exhibits  little  axial  variation  over  much  of  the 
cylinder  length  [ 3] ,  we  have  relaxed  the  boundary  conditions  at  the  end  walls. 

In  cylindrical  coordinates  r,9,z,  we  obtain  [3j  in  an  infinitely  long  cylinder 

V,  =  (0,0, watt.),  Pi  -  0  (2) 

where 

v,  “  ‘'.(f, 4)  **=  2<(/(r)eos^  +  j(r)sin*].  (3) 

and  g  and  /  are  the  real  and  imaginary  parts,  respectively,  of 

P(r)  ~  g  +  if  =  r  -  /,( qr)/I \(q).  (4) 

where  I\  is  the  modified  Bessel  function,  and  q  —  (1  4-  i)  (Pe/2)1^5.  This  solution  is  valid  for  arbi¬ 
trary  Re  but  may  be  unstable  as  Re  exceeds  some  critical  value.  Comparison  of  the  velocity  distribution 
with  computational  results  (5]  for  Re  —  14.9  has  shown  excellent  agreement  in  the  center  section  of 
the  cylinder.  At  the  higher  Reynolds  number  Re  =  45.7,  the  agreement  is  still  satisfactory,  with  sys¬ 
tematic  deviations  due  to  a  weak  cellular  motion  not  yet  incorporated  in  our  analysis. 

3.  Momenta 

With  given  deviation  velocity  V<(  the  moments  on  a  finite-length  section  of  the  cylinder  can  be 
calculated.  We  consider  a  control  volume  R  (surface  5)  formed  by  the  solid  cylindrical  wall  and  liquid 
surfaces  at  both  ends.  Conservation  of  angular  momentum  requires 

M  +  £(rXF<)  =  ±  f  J  f(rxVi)pdR  +  ///  jrX(2n  xV^jpdR  (5) 

R  R 

+  /  /  (r  XV,)(V,-n)pdS+  /  /  (r  xV,)(V,-n)p</S 

s  s 


APPENDIX  A 


A-27 


where  n  is  the  outer  unit  normal.  On  the  left-hand  side,  M  is  the  resultant  torque  on  the  control 
volume.  The  second  term  accounts  for  the  moments  due  to  the  shear  force  F*.  This  term  vanishes 
along  the  solid  sidewall  and  the  contributions  from  both  ends  cancel.  On  the  right-hand  side,  the  first 
term  vanishes  for  steady  V,.  The  second  term  originates  from  Coriolis  forces  in  the  nutating  system. 
The  third  term  vanishes  since  V*  has  only  an  axial  component.  The  last  term  then  provides  the  net 
rate  of  angular  momentum  flux  through  the  control  surface. 

Substitution  of  V*  leads  to  the  following  expressions  for  the  components  of  M: 


AY,  — ■  m<(2n  «sin0)(wa)  m, ,  m,  =*  -  J  r*fdr 

0 

(6a) 

1 

AY,  —*  m,(2D  «sin0)(u>«)  m,,  m,  —  -  J  T3gdr 

0 

(6b) 

1 

AY,  ^(20  atinfi)3  r  .,  m,  =  f  r3Jdr  *=  -  m. 

(6c) 

0 


where  m<  is  the  liquid  mass  in  the  cylinder.  In  this  form,  the  components  AY,,  AY,  represent  the  net 
rate  of  angular  momentum  flux  through  the  liquid  endwalls,  whereas  the  roll  moment  AY,  is  solely  due 
to  Coriolis  forces. 

A  different  interpretation  can  be  derived  using  the  differential  equation  for  F  [4]: 

m.  =  -  m,  -  /  r1  Jdr  =  -  »»,“-/  ^gdr  -  ^ ^  y-  (7) 

In  this  form,  the  moments  are  directly  related  to  the  shear  forces  at  the  cylindrical  sidewall,  r  *=  1. 
Since  /'(1)<0,  j'(l)<0,  the  roll  moment  AY,  is  always  positive  (even  for  D  <0).  From  (4),  we 
obtain 


F'(D  -r(l)+  *7  '( 1)  -  2  -  */. o(9 )//,(?)•  (8) 

For  small  Reynolds  numbers,  use  of  series  expansions  for  the  modified  Bessel  functions  provides  the 
approximations  m,  f=»/?e/98,  m,  -  /?ea/l!>36,  that  can  be  used  for  quick  estimates  up  to  Re  <10. 
These  results  disregard  the  effect  of  the  solid  end  walls  of  the  cylinder,  where  the  axial  flow  reverses 
direction.  We  expect  that  the  effect  of  this  flow  reversal  will  be  primarily  on  the  pitch  moment  AY,.  We 
also  expect  that  the  result  for  the  roll  moment  AY,  is  rather  accurate,  and  slightly  overestimates  the 
effect  of  Coriolis  forces.  Experimental  data  for  a  verification  of  our  results  are  scarce.  Only  the  roll 
moment  has  been  reliably  measured  in  spin-down  experiments  with  a  full-scale  cylinder  (lj. 

41.  Comparison  of  the  Roll  Moment  with  Experiments 

Figure  2  shows  the  comparison  of  theoretical  and  experimental  results  for  the  roll  coefficient  m, 
on  a  doubly  logarithmic  scale.  The  initial  spin  rate  w  =  4000  rpm  has  been  used  for  reducing  the 
experimental  data,  for  He  <  10,  the  experimental  data  match  the  analytical  data  ss  well  as  the  asymp¬ 
totic  law  m,  *«/?«/ 96.  The  coefficient  m,  assumes  a  pronounced  maximum  at  Re  19.  Whereas  this 
maximum  was  earlier  thought  to  originate  from  hydrodynamic  instability  with  respect  to  a  cellular 
motion,  we  find  a  simple  explanation  in  the  properties  of  the  axial  velocity  component  /  in  the 
1  ,t~  plane.  The  systematic  deviation  for  Reynolds  numbers  20  <  Re  <  200  may  be  attributed  to  either 
the  effect  of  a  cellular  motion  or  to  the  neglect  of  thr  end  walls.  The  two  data  points  at  Re  >  10s  are 
likely  to  be  for  a  turbulent  internal  flow. 
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Figure  1.  Definition  sketch 


Figure  2.  Comparison  of  the  theoretical  result 
for  m,  vs.  Re  with  experimental  data  ( 1 ) .  The 
stxajght  line  shows  the  asymptotic  law 
m,  ~  Re  /06. 


Figure  3.  Comparison  of  the  theoretical  result 
for  the  roll  moment  M,  vs.  kinematic  viscosity 
v  with  experimental  data  ( 1 ) .  Parameters: 
a  «=  60.3mm,  c/a  «=  4.29,  9  <=  20°, 

fl  ”  500rpm. 
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Figure  4.  Theoretical  results  for  the  roll  Figure  5.  Theoretical  results  for  the  roll 

moment  M,  vs.  kinematic  viscosity  u  for  moment  M,  vs.  nutation  rate  0  foi  different 

different  spin  rates  u .  Parameters  as  in  Fig.  3.  spin  rates  w.  Parameters:  a  =  50.4mm, 

e/a  —  4.5,  8  =  20°,  u  =  10* eS. 


Figure  6.  Theoretical  results  for  the  roll  Figure  7.  Experimental  results  for  the  roll 

moment  M,  vs.  spin  rate  w  for  different  moment  M,  vs.  spin  rate  cj  for  different 

kinematic  viscosities  v.  Parameters:  kinematic  viscosities  v.  Parameters: 

o  —  50.4mm,  c/a  *=  4.5,  8  =  20°,  a  =*■  50.4mm,  c/a  =  4.5,  *  —  20°, 

A  -=  625rpm.  O  «  600rpm. 
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ABSTRACT 


Spin-stabilized  projectiles  with  liquid  payloads  can  experience  a  severe  flight  instability  charac¬ 
terized  by  a  rapid  yaw-angle  growth  and  a  simultaneous  loss  in  spin  rate.  Laboratory  experi¬ 
ments  and  field  tests  have  shown  that  this  instability  originates  from  the  internal  fluid  motion  in 
the  rang  f  high  viscosity.  After  evaluation  of  the  experimental  data  and  analysis  of  the  equa¬ 
tions  for  it  fluid  motion  in  a  spinning  and  nutating  cylinder,  we  have  developed  a  simple 
model  of  this  flow.  Disregarding  the  finite  length  of  the  cylinder,  this  model  provides  the  flow 
field  and  the  viscous  contribution  to  the  liquid  momenta  in  analytical  form.  At  low  Reynolds 
number,  the  flow  field  sgrees  well  with  computational  results  for  the  center  section  of  a  cylinder 
of  aspect  ratio  4.3.  Ihe  roll  moment  caused  by  this  flow  largely  agrees  with  experimental  data 
for  a  wide  range  of  Reynolds  numbers.  Estimates  of  the  temperature  variation  indicate  that 
discrepancies  at  very  low  Reynolds  numbers  may  originate  from  associated  changes  of  the 
viscosity  during  the  experiments. 
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1.  Introduction 

It  is  well-known  that  spin-stabilised  shells  carrying  liquid  payloads  can  suffer  dynamical 
instability.  For  cylindrical  cavities  and  low  viscosity  of  the  liquid,  the  instability  due  to  basically 
inviscid  inertial  waves  can  be  predicted  by  the  Stewartson-Wedemeyer  theory  (Stewartson  1959; 
Wedemeyer  1959).  This  theory  rests  on  the  boundary-layer  approach  and  is,  therefore,  res¬ 
tricted  to  the  range  of  sufficiently  large  Reynolds  numbers.  The  instability  of  certain  shells  like 
the  XM761  (D’Amico  1977;  1978),  however,  escapes  such  a  prediction  and  is  also  dis¬ 
tinguished  in  character  owing  to  the  rapid  loss  in  spin  rate.  Experiments  with  a  full-scale 
liquid-filled  cylinder  (Miller  1982)  and  subsequent  field  tests  (D'Amico  Sc  Miller  1979)  estab¬ 
lish  that  this  new  flight  instability  is  most  pronounced  for  liquid  fills  of  very  high  viscosity. 

We  conduct  a  theoretical  analysis  of  this  problem  in  order  to  support  the  ongoing  experi¬ 
ments  and  to  independently  obtain  insight  into  the  anatomy  of  the  flow  phenomena.  The  initial 
steps  of  this  analysis  are  reported  elsewhere  (Herbert  1982):  evaluation  of  the  experimental 
data  base,  dimensional  analysis,  scaling  aspects,  governing  equations,  and  discussion  of  various 
simplifying  assumptions.  Two  observations  in  this  earlier  work  led  to  the  approach  discussed  in 
the  following.  First,  if  the  despin  (negative  roll)  moments  (Miller  1982)  and  void  observations 
(Miller  1981)  are  correlated  with  the  Reynolds  number  Re,  at  least  three  regions  can  be  dis¬ 
tinguished.  At  low  Re,  the  despin  moment  increases  proportional  to  Re,  and  the  void  in  an 
incompletely  filled  cylinder  is  parallel  to  the  spin  axis.  This  suggests  a  simple  fluid  motion  that 
is  essentially  independent  of  the  axial  coordinate,  except  in  the  neighborhood  of  the  end  walls. 
In  a  middle  range  of  Re ,  the  despin  moment  assumes  a  maximum,  and  a  wavy  distortion  of  the 
void  seems  to  indicate  a  cellular  structure  of  the  fluid  motion.  This  cellular  motion  can,  in 
principle,  originate  from  hydrodynamic  instability  of  the  basic  flow  with  respect  to  axially 
periodic  disturbances.  At  still  higher  Reynolds  numbers,  the  despin  moment  decreases  «ith 
increasing  Re  in  a  manner  not  clearly  defined  by  the  few  available  data  points.  The  void  obser¬ 
vations  indicate,  however,  that  the  motion  ultimately  becomes  turbulent. 

Tie  second  observation  is  the  appearance  of  the  nutation  rate  and  angle  as  a  small  param¬ 
eter  in  the  equations  for  the  deviation  from  solid-body  rotation.  The  forcing  term  due  to  nuta¬ 
tion  can  be  considered  small  enough  for  linearization  of  the  equations  in  the  situations  of  prac¬ 
tical  interest. 

In  the  following,  wc  describe  the  development  of  a  simple  system  of  equations  for  the 
basic  flow.  Analytical  solutions  are  given  for  the  flow  field,  for  the  liquid  moments,  and  for  the 
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rate  of  change  of  temperature.  A  comparison  is  made  with  computer  simulations  of  the  flow 
(Vaughn  et  al.  1983;  198S)  and  with  experimental  data  for  the  moments  (Miller  1982). 

2.  Governing  Equations 

We  consider  the  motion  of  a  fluid  of  density  p  and  viscosity  p  in  a  cylinder  of  radius  a 
and  length  2c  that  rotates  with  the  spin  rate  u  about  its  axis  of  symmetry,  the  z-axis.  We  con¬ 
sider  the  motion  with  respect  to  the  nutating  coordinate  system  *,y,r.  This  system  is  obtained 
from  the  inertial  system  X,Y,Z,  by  a  rotation  with  the  nutation  angle  0  about  the  axis  Y=y. 
Therefore,  z  is  in  the  Z,»- plane,  and  this  plane  rotates  about  the  J-axis  with  the  nutation  rate 
Q  .  The  two  axes  of  rotation  intersect  in  the  center  of  mass  of  the  cylinder,  as  shown  in  figure 
1.  In  contrast  to  the  experimental  procedures  (Miller  1C82),  we  consider  w>0,  fl,  and 
0<6<ir /2  as  constant. 

The  fluid  motion  is  governed  by  the  Navier-Stokes  equations  written  in  the  nutating  coor¬ 
dinate  system: 

+  20  X  V. +  n  X  (flXr)]  -  -?P.+fiV*V.  ,  (la) 

v  v„  *=  0  .  (lb) 


V.  is  the  velocity  measured  in  the  nutating  frame,  Pm  the  pressure,  and  r  the  position  vector. 
The  body  force  due  to  gravity  has  been  disregarded.  Equations  (1)  are  subject  to  the  no-slip 
and  no-penetraUon  conditions  at  the  cylinder  walls. 

It  is  convenient  (Herbert  1982)  to  split  the  velocity  and  pressure  fields  according  to 

v.  =  v,  +  V,  ,  P,  =  P,  +  Pt,  (2) 


where  Y,,P,  describe  the  state  of  pure  solid  body  rotation,  whereas  represent  the  devia¬ 

tion  from  solid  body  rotation.  The  advantage  of  this  isolated  view  on  the  deviation  is  obvious: 
V,j  and  the  reduceJ  pressure  Pi  are  responsible  for  the  observed  flight  instability.  A  glance  at 
the  equations  shows  that  V4  =  0  and  Pd  e  0  if  either  one  of  the  following  conditions  is 
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satisfied:  w  =  0,  fi  *=  0,  0  =  0  or  p  -»  oe  (solid  fill). 

The  equations  for  N4,P4  are  written  in  terms  of  nondimensional  quantities  v* ,p4.  We  use 
a,ui,  and  p  fo*  sealing  length,  time  and  mass.  Note  that  this  choice  is  ambiguous  (Herbert 
1982)  and  excludes  the  case  w=0  which  lacks  practical  interest.  The  problem  then  depends  on 
four  nondimensional  parameters: 

>  =  e/a  aspect  ratio 

6  nutation  angle 

t  =  fi  /w  frequency 

Re  =  puai /(i  Reynolds  number. 

The  aspect  ratio  enters  the  solution  only  through  the  boundary  conditions.  The  boundary  con¬ 
ditions  on  vj  are  homogeneous. 

In  cylindrical  coordinates  r,^,z,  the  equations  for  the  nondimensional  deviation  velocity 
vs  —  (tV.ti^u,)  and  pressure  p4  take  the  form 


Id.  1  dv4  dv, 

r  dr  n,r  r  d<t>  dz  ”  ' 


D'vr  - - 2(1  +  t,)v4  +  2 t4v,  = 


r2  r2  d*  1  ’ 


V.Vi 

D  v4  +  -  +  2(  1  +  r,)vr  -  2 rrv,  = 


1  dpt  i 


r  d<t>  +  Rc^D"V*  ~  -2  +  -2 


2  dvr. 


r*  r*  dd> 


where 


dpj  1 

D'v,  +  2rrv4  -  2r4vr  =  -  -  2rrr  +  — £>%,  . 


m  d  ^  d  ^  d  ^  d 

D  =  Tt  +  1*  ~r  +  w*  57  ’ 


Z>"  = 


a2  i  a  i  a2  a2 

ar2  +  r  dr  ^  r2  a*2  +  ar2  ’ 


(3a) 

(3b) 


(3c) 


(3d) 


APPENDIX  A 


Htrhert 


4 


A~ 


and 


r,  =  -ecos^  ,  r4  =*  esind  ,  r,  *  rcosd  ,  t  ■»  rsintf  . 


(<) 


The  primary  effect  of  nutation  is  contained  in  the  ^-periodic  iorce  term  -2 rrr  •»  2<rcosd  in 
the  /-momentum  equation  (3d).  If  this  term  vanishes  throughout,  t  “  0,  equations  (3)  sup¬ 
port  a  trivial  solution  vd  m  0,  pt  m  0. 

The  system  (3)  of  equations  is  similar  to  the  system  numerically  solved  by  Vaughn  et  al. 
(1983;  1983),  but  simplified  by  introducing  the  reduced  pressure  pd.  We  also  note  that  this  sys- 
tem  supports  certain  symmetries.  Let  v„v4,v,  and  pd  be  the  solution  at  point  r,d,r,  then  the 
velocities  and  pressure  at  the  corresponding  point  r,d  +  w,-  z  are  v,,v4l-  v,  and  p4.  These  sym¬ 
metries  can  be  exploited  for  essential  savings  in  computational  work. 

2.1  Linearized  equation* 

For  sufficiently  small  t  yd  0,  it  is  obvious  that  the  deviation  velocity  is  of  order  0(e).  In  the 
situations  of  practical  interest,  <  °*=  (0  /w)sin0  turns  out  to  be  a  rather  small  parameter.  Even 
a  conservative  estimate  with  O  <  500  rpm,  w>  3000  rpm,  and  $<  20*  provides  values  of 
«<  0.054.  Consequently,  it  seems  well  justified  to  linearise  the  equations  in  e.  This  linearisa¬ 
tion  imposes  no  restriction  on  the  Reynolds  number. 

While  the  continuity  equation  remains  unaffected,  linearisation  of  the  momentum  equa¬ 
tions  provides 


D'vr  -  2(1  +  r,)v,  = 


D'v4  +  2(1  +  r,)vf  — 


D'v,  =  - 


_ dp*  +  1  {D"V  _  .  2  dv*> 
dr  Re[  '  r3  r2  d*  1  ’ 

(5a) 

'//; +  s  +  jiji- 

(5b) 

+  -kD%‘  ■ 

(5c) 

The  system  (3a),  (6a)-(5c)  of  equations  is  still  quite  difficult  to  solve.  Any  serious  attempt  to 
s&usfy  all  boundary  conditions  leads  directly  to  a  purely  computational  approach.  Use  of  the 
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boundary-layer  approximation  would  aimplify  the  task  but  seems  inappropriate  in  the  interesting 
range  of  low  Reynolds  numbers. 
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3.  The  core  flow 

We  recall  that  the  flow  in  a  relatively  long  cylinder  (aspect  ratio  X  =  4.3)  at  low  Reynolds 
number  is  expected  to  have  a  rather  simple  structure  and  to  provide  a  roll  moment  proportional 
to  Re  (Herbert  1982).  Closer  analysis  of  the  equations  suggests  that  this  flow  exhibits  little 
axial  variation  over  much  of  the  cylinder  length.  The  effect  of  the  end  walls  will  be  essential 
only  over  an  axial  distance  of  0(1)  from  the  ends.  Therefore,  we  have  relaxed  the  boundary 
conditions  at  the  end  walls.  In  this  way,  we  seek  a  steady  flow  in  a  finite  segment  of  an 
infinitely  long  cylinder. 

The  r-independent  force  term  in  eq.  (5c)  can  be  balanced  only  by  a  purely  axial  deviation 
velocity.  It  is  consistent  with  the  linearized  equations  to  assume  a  solution  in  the  form 

v,  =  (0,0, v,)  ,  Pi  =  0  .  (6) 

Moreover,  since  v,  is  of  order  0(<)  and  periodic  in  d>,  we  write 

t>,  =  v,(r,<t>)  =  2<(/(r)cos^  +  f(r)sin^)  ,  (7) 

where  /  and  g  are  the  imaginary  and  real  parts,  respectively,  of  the  complex  function 

f(  *■)  =  «('’)  +  '/(*■)  (8) 

Substituting  (6)-(8)  into  the  linearized  equations  and  the  no-slip  conditions  at  the  cylinder  wall 
provides 


r2F"  +  r  F'  -  (1  +  iRt  r 3)F  = 
F  =  0  at  r  =  1  , 

F  finite  at  r  =  0  , 


i  Re  r3  ,  (9a) 

(9b) 

(9c) 


where  (9c)  is  necessary  for  a  physical  solution.  The  primes  denote  d/dr. 


3  1  Solution  for  Re  — < >0  and  Rt—o o 
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For  Re~*Q,  the  solution  of  equations  (8)  equations  can  be  found  in  the  form  of  series 
expansions  in  Re , 


/  =  x(f-rS>-  ?5T?(7r"  12S  +  6r‘-  r>)  +  0(Rc>) , 


9  =  ik[2r  _  3fS  +  ri)  +  °(Rei) 


(10a) 

(10b) 


With  higher  terras  included,  these  series  converge  for  Re<  12. 

In  the  limit  Re-* oo,  one  obtains 

/— *0,  g  —•  r  as  Re  —•  oa.  (11) 

Owing  to  the  loss  of  the  highest  derivatives,  however,  this  solution  cannot  satisfy  the  boundary 
conditions  (8b)  and  is  valid  only  outside  thin  boundary  layers  near  the  wall  at  r  =  1. 

Even  without  any  knowledge  of  the  solution  in  the  intermediate  range,  the  different  char¬ 
acter  of  the  basic  flow  at  low  and  high  Reynolds  numbers  is  evident.  At  low  Re ,  the  component 
/  in  the  *,*-plane  <f>  =  0  dominates  the  solution.  At  high  Re,  f  is  negligible  except  near  the 
wall  of  the  cylinder  while  g  in  the  y,z-plane  d>  —  80*  is  dominating.  One  might  well  expect 
that  the  initial  linear  increase  of  /  with  Re  and  the  change  in  the  flow  structure  is  related  to  the 
observed  properties  of  the  roll  moment. 

3.2  Solution  for  arbitrary  valuta  of  Re 

In  earlier  work  (Herbert  1983),  we  have  applied  a  spectral  collocation  method  for  numeri¬ 
cally  solving  a  real  system  of  equations  for  /  and  g  equivalent  to  eqs.  (9).  Series  in  odd  Che- 
byshev  polynomials  for  the  interval  0  <  r  <  1  provide  accurate  solutions  at  rather  low  trunca¬ 
tion.  This  experience  together  with  the  minor  effect  of  harmonics  in  the  azimuthal  direction  at 
small  (  suggests  the  use  of  spectral  methods  for  efficiently  solving  the  nonlinear  equations  (3). 

Here,  we  derive  an  analytical  solution  for  the  core  flow  in  a  sufficiently  long  cylinder.  A 
particular  solution  of  the  inhomogeneous  equation  (9a)  is  F0  =  r,  whereas  the  homogeneous 
part  of  (9a)  is  the  equation  for  the  modified  Bessel  functions  /f(9r)  and  Ki(7r)  of  the  com¬ 
plex  argument  gr  where  g  =  (1  +  i)(  Re  /2)1''3.  In  order  to  satisfy  (9c),  K,(gr)  cannot  contri¬ 
bute  to  the  solution.  Finally,  (9b)  provides 
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F(r)  =  »  +  •'/  =  r  -  /i(?r)/A(?)  • 


(12) 


This  solution  is  valid  for  arbitrary  Rt  but  may  be  unstable  as  Re  exceeds  some  critical  value. 
Although  expressible  in  simple  form,  the  resulting  flow  field  exhibits  very  interesting  proper¬ 
ties. 

Rewriting  the  solution  in  terms  of  Kelvin  functions  of  real  argument  is  of  little  advantage 
for  the  numerical  evaluation.  We  have  used  a  combination  of  ascending  series  and  asymptotic 
expansions  for  large  arguments  (Abramowiti  &  Stegun  1672)  for  evaluating  F(r).  With  the 
solution  (12)  at  hand,  it  is  straightforward  to  derive  the  approximations  (10)  from  the  ascend¬ 
ing  series  for  /,  (and  to  explain  the  convergence  problem  for  larger  Re).  Complementary  to 
(11),  the  asymptotic  expansion  for  large  arguments,  i.e.  large  Reynolds  numbers  provides  the 
boundary-layer  behavior 


f«r-  . 


(13) 


This  expression  agrees  to  within  l%with  (11)  provided  r<l-$.  The  boundary  layer  thickness 
6  can  be  obtiined  from  the  transcendental  equation 


S  =  \flJWe  (4.605  -  -|ln(l  -  «))  , 


(H) 


e.g.,  6  =  0.223  for  Re  =  1000. 


3,3  77i«  velocity  field 

We  have  chosen  three  different  graphical  representations  in  order  to  illustrate  the  charac¬ 
teristic  changes  of  the  velocity  distribution  over  the  cylindrical  cross  section  with  increasing  Re . 
Figure  2  shows  the  components  /  (in  the  r,r-plane)  and  g  (in  the  y,r-plane)  for  a  wide  range 
of  Reynolds  numbers.  The  opposite  sign  of  the  velocity  at  diametral  points  assures  xero  net 
flux  of  mass  through  the  cross  section.  The  curves  represent  cuts  through  the  contour  plots  of 
these  functions  of  r  and  Re  in  figure  3  at  the  tick  marks  Re  =1,  10, 100,and  1000.  Up  to 
Rc  ~  5,  the  velocity  distribution  is  governed  by  /.  This  component  n<>v«.r  fVffeH.i  a  value  of 
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0.4,  usumes  a  maximum  at  i?e  *=»  20  and  retains  significant  site  only  in  a  shrinking  neighbor¬ 
hood  of  the  wall  as  Re  increases.  The  component  g  rapidly  increases  from  negligible  values  as 
Pe  >  S  and  approaches  the  linear  increase  with  r  according  to  (11)  except  near  the  wall  at 
r  ==  I.  In  figure  4,  the  data  of  figure  2  are  combined  into  contour  plots  of  the  axial  velocity 
vJ/(2e)  over  the  cylindrical  cross  section.  These  plots  clearly  show  the  shift  of  the  velocity  max¬ 
imum  (marked  by  -f )  from  ^~0at/?e  =  ltod’  ~C0*  at  Re  =  1000.  Figure  4d  also  illus¬ 
trates  the  ramp-like  velocity  distribution  over  most  of  the  cross  section  and  the  boundary  layers 
with  6  «=  0.223. 

Superposition  of  the  deviation  velocity  Vg  and  the  solid  body  rotation  V,  according  to  eq. 
(2)  leads  to  an  asimuthally  periodic  velocity  field  V,  which  is  steady  in  the  nutating  frame.  The 
paths  of  fluid  elements  are  circular  orbits  about  axes  that  ue  inclined  to  the  ;-axis.  The  incli¬ 
nation  depends  on  radius  and  Reynolds  number. 

Figure  5  compares  the  dimensional  velocity  distributions  obtained  from  (7),  (12)  with 
computational  results  for  the  center  cross-section  ( t  =  0)  of  a  cylinder  of  aspect  ratio  4.3.*  The 
agreement  for  Re  ■=  14.9  is  considered  representative  for  the  range  of  lower  Reynolds 
numbers.  We  have  repeated  the  numerical  simulation  of  the  Sow  at  this  Reynolds  number  with 
a  modified  version  of  the  Saadia  code  and  obtained  very  small  components  |vr  |  <  0.005  m/s, 
|tig  |  <  0.05  m/s  at  *  =  0.  These  results  verify  our  estimates  and  justify  the  use  of  linearised 
equations.  Moreover,  disregarding  the  presence  of  end  walls  seems  to  have  little  effect  in  the 
center  portion  of  the  cylinder.  The  radial  distribution  of  V,  in  the  range  -3.5  <  *  <3.5  is 
nearly  identical  with  the  data  shown  in  figure  5. 

Figure  8  shows  a  similar  comparison  for  Re  =  45.7.  At  this  higher  Reynolds  number,  we 
find  a  systematic  deviation  between  the  theoretical  result  and  numerical  results  at  different  axial 
positions.  We  attribute  this  deviation  to  a  superposed  cellular  motion  that  is  not  yet  incor¬ 
porated  into  our  analysis. 


•The  data  were  kindly  provided  by  Dr.  H.  Vaughn,  Sandia  National  Laboratories. 


APPENDIX  A 


A-42 


Th.  Herbert 


10 


4.  Moments 

With  the  deviation  velocity  V4  =  (0,0, wav,)  and  v,  given,  the  momenta  on  a  finite- 
length  section  of  the  cylinder  can  be  calculated.  We  consider  a  control  volume  R  (surface  5) 
formed  by  the  solid  cylindrical  wall  and  liquid  surfaces  at  both  ends.  Conservation  of  angular 
momentum  requires 


M  +  £(rxF,)=  ±J  /  f(rxVd)pdR  +  /  /  /  |r  X(20  xV4)]ptt  (15) 

R  R 

+  //(rXV,)(V,  n)pdS+//(rXV,)(V<  n)pdS  , 
s  s 


where  n  is  the  outer  unit  normal.  On  the  left-hand  side,  M  is  the  resultant  torque  on  the  con¬ 
trol  volume.  The  second  term  accounts  for  the  moments  due  to  the  shear  force  and  van¬ 
ishes  owing  to  the  solid  sidewall  and  cancellation  of  the  contributions  from  both  ends.  On  the 
rightrhand  side,  the  first  term  vanishes  for  steady  V*.  The  second  term  originates  from 
Coriolis  forces  in  the  nutating  system.  The  third  term  vanishes  since  V,  has  only  an  axial  com¬ 
ponent.  The  last  term  then  provides  the  net  rate  of  angular  momentum  flux  through  the  con¬ 
trol  surface. 

Substitution  of  V*  leads  to  the  following  expressions  for  the  cartesian  components  of  M: 


M,  =  m<(2n  a  sin0)(  wa )  m,  ,  m,  —  -  J  r“fdr  ,  ( 16a) 

o 

i 

A/r  =  m, { 20  asinfl)(wa)  mf  ,  mf  =  -  J  r2gdr  ,  ( 16b) 

o 


A/,  —  mj(20  aaintf)2  m,  ,  m,  —  J  r2  J  dr  =  -  m,  ,  (16c) 

0 


where  m*  =  2npa7e  is  the  liquid  mass  in  the  cylinder.  In  this  form,  the  components  A/,,  Mt 
represent  the  net  rate  of  angular  momentum  flux  through  the  liquid  ends,  whereas  the  roll 
moment  M,  is  solely  due  to  Coriolis  forces.  A  close  relation  between  roll  moment  A/,  and  yaw- 
moment  SI,  has  also  been  found  by  Murphy  (1984,  1685) 
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A  different  interpretation  can  be  derived  using  the  differential  equation  (9a),  integrating 
by  parte,  applying  (9b),  and  separating  real  and  imaginary  parts: 


f  2  j  /'(l)  1 

rtdr  =  -Rr-  7 


(17a) 

(17b) 


In  this  form,  the  momenta  are  directly  related  to  the  shear  forces  at  the  cylindrical  sidewall, 
r  =  1 .  Since  /  '( 1)  <  0,  g  '( 1) <  0,  the  roll  moment  M,  is  always  positive  (even  for  0  <  0), 
while  M,  is  negative  for  fi  >  0  and  changes  sign  with  fl  .  For  small  Re,  the  series  (10)  pro¬ 
vide  the  approximations 


m. 


Rt2 
1536  ’ 


( 18) 


that  can  be  used  for  quick  estimates  up  to  Re<  10.  The  linear  increase  of  m,  and  M,  with  Re  is 
consistent  with  the  experimental  data.  From  the  analytical  solution  (10),  we  obtain 


F'(  1)  =  f'(l)  +  ./'(l)  =  2-  9/oU)//»U)  (I0) 

Substitution  into  (17)  provides  the  variation  of  m,,  m,  with  the  Reynolds  number  shown  in 
figure  7.  The  coefficient  m,  a-  >rnes  a  pronounced  maximum  at  Re~l9.  The  occurrence  of 
this  maximum  was  earlier  thought  to  originate  from  hydrodynamic  instability  with  respect  to  a 
cellular  motion.  Here,  we  find  a  simple  explanation  in  the  properties  of  the  axial  velocity  com- 
,.oiicnt  /  m  the  z,z-  plane  and  the  derivative  p'(l)  The  coefficient  ni^  is  negligible  for 
Re  <  5,  '  vrply  decreases  with  increasing  fie  ami  reaches  an  asymptotic  value  of  mf  -  1/4 
as  Re  -*  oo.  Hem  c,  for  (1  >  0,  A/f  tends  to  reduce  the  pitch  moment  due  to  the  solid  body 
if,  ition.  We  note,  however,  that  these  moments  represent  only  the  effect  of  viscous  shear  at 
the  cylindrical  side  wall  Shear  at  the  end  walls  and  the  contribution  of  the  pressure  are 
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neglected. 

The  data  base  for  the  yaw  and  pitch  momenta  is  scarce.  Computations  by  Vaughn  et  al. 
(1985)  indicate,  however,  that  the  pressure  contributions  to  these  momenta  are  larger  (and 
opposite  in  sign)  than  the  viscous  components.  Only  the  viscous  component  can  be  estimated 
from  our  solution.  Therefore,  we  concentrate  in  the  folio*  'ng  on  *  detailed  comparison  for  the 
roll  moment. 

In  figure  8  we  compare  the  asymptotic  law  (18)  and  the  theoretical  result  (17)  with  exper¬ 
imental  data  (Miller  1682)  and  computational  results  (Vaughn  et  al.  1685)  for  the  roll 
coefficient  m,  on  a  doubly  logarithmic  scale.  The  initial  spin  rate  u  >»  4000  rpm  has  been  used 
for  obtaining  the  nondime naional  values  from  the  experiment.  For  Re  <  10,  the  experimental 
data  match  the  analytical  result  as  well  as  the  asymptotic  law  m,  hie/98.  The  deviation 
between  theoretical  and  computational  results  is  probably  due  to  a  larger  axial  extend  of  the 
end  effects  at  very  low  Reynolds  numbers.  Good  agreement  with  the  computational  results  is 
obtained  near  the  maximum  of  m,.  The  point  at  Re  ■—  113  is  close  to  the  Reynolds  number 
where  the  numerical  simulation  fsils  to  converge  to  a  steady  solution,  and  may  not  be  very 
accurate  The  experiments  find  the  maximum  roll  moment  at  (lightly  lower  Reynolds  numbers 
than  the  theoretical  value.  Iu  fact,  this  discrepancy  will  increase  as  lower  spin  rates  w  are  used 
for  data  reduction  In  view  of  the  agreement  between  theoretical  and  computational  results,  the 
discrepancy  can  not  arise  from  the  approximations  employed  in  our  analysis.  A  first  possible 
source  may  be  the  effect  of  unsteadiness  in  the  spin-down  experiments.  More  likely,  however, 
the  shift  is  can-  i  by  changes  of  temperature  and  viscosity  during  the  experiments.  A  moderate 
increase  in  temperature  would  reduce  the  viscosity  of  the  working  fluids  (silicone  oii,  corn 
syrup)  and  hence  shift  the  maximum  to  higher  Reynolds  numbers.  Miller  (personal  communi¬ 
cation)  observed  a  temperature  increase  by  as  2.6'C  per  run  up  to  ss  10*C  above  ambient  tem¬ 
perature  after  repeated  runs.  Vaughn  et  al.  (  1685)  used  these  values  for  correcting  the  result*, 
witli  some  improved  agreement.  We  waive  such  a  correction  but  discuss  the  temperature 
increase  in  more  detail  in  the  next  chapter. 

Am  a  final  observation  in  figure  8,  we  note  the  change  in  tendency  for  the  twc  experimen¬ 
tal  data  points  at  R<  >  105.  It  is  likely  that  the  internal  flow  become,  unsteady  and  ultimately 
turbulent  as  the  Reynolds  number  increases.  Preliminary  result*  from  flow  visualiiation  in  * 
sniall-icaJe  experiment  (Pierpond  1085)  indicate  that  these  two  point*  are  for  a  turbulent  inter¬ 
nal  fluv, 
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In  figure  9  we  recast  experimental,  computational  and  theoretical  results  for  the  dimen¬ 
sional  roll  moment  M,  in  different  form.  Whereas  the  asymptotic  properties  are  concealed,  the 
linear  scale  for  M,  reveals  the  pronounced  maximum  of  the  roll  moment  for  viscosities  near 
v  =  106  cSl  and  more  clearly  indicates  that  theory  and  computation  yield  larger  maximum 
values  than  the  despin  experiments  with  the  old  test  fixture  (Miller  1982).  More  recent  meas¬ 
urements  with  a  new  test  fixture  at  higher  spin  rates  (Miller,  personal  communication)  provide 
larger  maximum  values  slightly  in  excess  of  the  theoretical  result. 

For  the  roll  moment  as  a  function  of  nutation  angle  and  rate,  Herbert  (1983)  derived 
from  Miller’s  data  (1982,  fig.  12)  the  empirical  relation  M,  =  0.00814  (Osin#)2  Nm.  The 
theory  provides  M,  in  the  same  form  but  with  a  somewhat  larger  factor  of  0.0111.  This  com¬ 
parison  for  a  fluid  of  kinematic  viscosity  v  =»  210s  cSt  is  likely  to  be  biased  by  temperatuie 
effects.  A  notable  feature  of  the  roll  moment  as  a  function  of  nutation  rate  at  different  spin 
rates  is  shown  in  figure  10.  For  these  parameters  in  the  range  of  the  maximum  roll  moment, 
the  dependence  of  M,  on  u  is  non-monotonic,  e  g.  the  data  for  lj  =*  9000  rpm  are  in  between 
those  for  w  =  3000  and  9000  rpm.  This  puzzling  behavior  has  been  observed  by  Miller  in 
experiments  with  the  new  test  fixture.  From  the  theoretical  result  it  is  obvious  that  M, 
decreases  (increases)  with  w  for  sufficiently  small  (large)  viscosities  to  the  left  (right)  of  the 
maximum  in  figure  9. 

The  interpretation  of  the  experimental  results  has  been  hampered  by  the  observation  of 
Miller  (1982)  thui  “the  despin  moment  was  not  s  function  of  the  canister  spin  rate,  provided  a 
sufficient  spin  rate  w  present'-.  In  contrast,  the  theoretical  result  (14c),(15a),(19)  depends  on 
the  spin  rate  since  q  ~  Re 1/15  and  Re  —  w  for  fixed  a  and  w.  Figure  11  shows  the  tneoretical 
results  for  M,  as  a  function  of  the  spin  rate  ui  for  viscosities  v  «=  10s,  104,  and  10s  cSl  on  linear 
scales  Note  that  in  some  range  of  w,  M,  appears  indeed  nearly  independent  of  the  spin  rate, 
especially  for  u  =  106  cSt  where  the  maximum  of  A/,  stretches  out  over  most  of  the  observed 
range  (3000  <  w  <  9000  rpm)  of  spin  rates.  Figure  11  also  shows  different  prototypes  of 
behsvor  'hat  &r>,  distinguished  by  the  position  of  the  maximum  roll  moment  along  the  ui  axis 
Experimental  data  for  similar  conditions  are  shown  in  figure  12  and  verify  the  theoretically 
predicted  behavior.  Moreover,  these  data  suggest  major  simplifications  in  the  experimental  pro¬ 
cedures  Whereas  the  experimental  data  in  figuie  9  were  obtained  by  using  numerous  working 
fluids  of  different  viscosities,  a  more  complete  set  of  data  can  be  generated  by  carefully  moni¬ 
toring  the  spin-down  for  s  few  runs  with  fluids  in  the  range  of  low,  medium,  and  high  viscosi¬ 
ties  as  in  figure  12 
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S.  Temperature  effect 

The  comparison  of  theoretical  and  computational  results  with  experimental  data  seems  to 
be  biased  by  the  effect  of  increasing  temperature  on  the  viscosity  of  the  working  fluid.  These 
effects  appear  more  pronounced  at  high  viscosities  and  high  spin  rates.  For  an  estimate  of  the 
rate  of  change  of  the  average  temperature  T,  we  consider  a  control  volume  R  (surface  S) 
formed  by  liquid  surfaces  along  the  cylinder’s  side  and  end  walls.  The  material  properties  are 
assumed  to  be  constant  and  heat  transfer  through  the  surface  is  disregarded.  Balancing  the  rate 
of  change  of  energy  with  the  work  done  on  the  control  volume,  we  obtain  after  some 
simplifications 

IJr  V.dS,  (20) 


where  e,  is  the  specific  heat,  r  the  vector  of  tangential  stresses,  and  V,  the  velocity  measured 
in  an  inertial  frame.  Since  is  independent  of  i,  the  contributions  from  the  cylinder  ends 
cancel.  The  only  contribution  is  due  to  the  shear  stress 

r„  —  n  ^  1,-t  —  asin* |/'(l)cos*  +•  f'(l)sin<£]  (21) 

in  the  axial  direction.  The  relevant  axial  component  of  the  velocity  (U  +  w)  Xrof  some  point 
on  the  surface  S  is  given  by  -  ft  asinfain^.  Integration  over  the  cylindrical  surface  yields 

dT 

c,  -r-  —  -  p(20  a  sin0)Jrr  oe  g'(  I )  .  (22) 

at 


After  substituting  for  m(  and  introducing  the  Reynolds  number,  this  result  can  be  written  as 


dT 

dl 


=  — — ( 2fl  a sinS)2( 


g'(D 

Rt 


(23) 


Comparison  with  eq.  (15a)  shows  that  the  rate  of  change  oT  temperature  can  be  directly 
expressed  in  terms  of  the  roll  moment, 


£T 

dl 


2(1  a  sing)2  m, 

2c, 


w 

2  *»|C, 


M. 


(24) 


This  result  immediately  shows  that  tlir  temperature  rise  per  run  cannot  be  specified  as  a  single 
number,  nor  should  a  uniform  correction  be  applied  to  the  experimental  data  Moreover,  tbe 
temperature  changes  increase  with  the  spin  rate,  and  consequently  are  quite  different  l or  the 
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experiments  with  tie  old  (Miller  1982)  end  the  new  test  fixture.  Using  the  maximum  value 
m,  ss  0.0854,  we  obtain  for  the  1982  experiments  (a  —  60.3  Tnm,  w  =  4000  rpm,  f)  =»  500 
rpm,  9  =  20*)  with  silicone  oil  (e,  w  1600J/(kg*C)  a  temperature  rise  of 
dT/dt  =  0.055'C/s.  Using  the  same  fluid  in  the  new  test  fixture  (a  =  55.4  mm,  u  =  104 
rpm,  fl  —  600  rpm,  9  =  20')  leads  to  a  temperature  increase  of  dT/dt  —  0.158'C/s 

A  single  run  consists  of  three  phases  (Miller,  personal  communication).  The  spin-up 
period  of  as  30  s  is  followed  by  a  sudden  start  of  the  nutation  and  a  period  of  ~  30  s  in  order 
to  reach  steady  conditions.  Finally,  the  shutdown  of  the  spin  drive  is  followed  by  a  spin-down 
period  of  15  s.  The  second  period  at  nearly  steady  conditions  and  maximum  spin  rate 
appears  most  relevant  to  the  modification  of  the  viscosity.  At  the  start  of  the  third  phase,  the 
average  temperature  may  have  increased  by  ssl.56*C  in  the  1982  experiments,  and  by 
as  4.75' C  in  the  more  recent  experiments  at  higher  spin  and  nutation  rates  These  values  are 
for  conditions  of  maximum  roll  moment,  and  may  be  considerably  lower  in  other  cases.  The 
value  of  2.5'C  measured  in  the  new  fixture  is  well  within  the  estimated  range.  Our  result  (24) 
indicates,  however,  that  a  single  measurement  is  insufficient  for  evaluating  the  temperature 
effects,  especially  those  in  a  different  experimental  setup. 
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0.  Concluding  remarks 

We  have  developed  a  simple  model  of  the  viscous  fluid  motion  in  s  spinning  and  nutating 
cylinder.  The  disregard  of  the  end  walls  has  some  obvious  consequences:  the  turning  Sow  near 
the  ends  and  the  associated  contributions  of  pressure  and  shear  stresses  to  the  moments  cannot 
be  obtained  from  this  model.  Nevertheless,  we  gather  understanding  as  well  as  quantitative 
information.  The  velocity  field  of  the  core  flow  agrees  well  with  computational  results  for  low 
Reynolds  numbers.  The  analytical  result  is  an  evident  example  for  the  formation  of  boundary 
layers.  The  core  flow  can  be  utilized  as  a  basic  flow  in  studies  of  hydrodynamic  instability  with 
respect  to  cellular  motions.  The  parametric  excitation  of  such  cells  by  the  azimuthally  periodic 
deviation  has  been  discussed  by  Herbert  ( 1884} .  The  core  flow  also  represents  the  lowest-order 
approximation  to  the  solution  of  the  nonlinear  equations  (3)  and  can  be  extended  by  higher- 
order  terms  in  c , 

The  roll  moment  agrees  well  with  measured  and  computed  values,  and  can  also  be  found 
at  Reynolds  numbers  too  large  for  successful  numerical  simulations.  The  roll  moment  ori¬ 
ginates  from  Coriolis  forces.  While  the  direct  calculation  of  the  yaw  moment  Buffers  from 
neglecting  the  pressure  contribution,  the  yaw  moment  can  be  found  from  the  roll  moment 
using  the  relations  given  by  Murphy  (1984,  1985).  The  pitch  moment  remains  an  open  issue. 
The  estimates  for  the  change  in  average  temperature  need  further  verification  once  more 
detailed  data  become  available. 

The  simple  form  and  scaling  relations  of  our  results  provide  guidance  for  sorting  and 
evaluating  the  experimental  data  base.  The  results  also  suggest  various  improvements  in  the 
experimental  procedures.  First,  tne  changes  in  temperature  and  viscosity  should  be  carefully 
monitored.  With  the  effective  viscosity  known,  a  closer  agreement  between  theory  and  obser¬ 
vation  is  to  be  expected,  Second,  the  yet  neglected  variation  of  the  roll  moment  with  the  spin 
rate  is  considered  relevant  and  in  fact  provides  the  roll  moment  in  some  range  of  Reynolds 
numbers.  Instead  of  producing  the  data  for  figure  9  by  using  numerous  viscosities  at  fixed  spin 
rate,  very  similar  data  can  be  generated  by  varying  the  spin  rate  for  a  few  fluids. 
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FIGURE  CAPTIONS 


Definition  sketch 

Components  /  and  g  of  the  axial  velocity  v,/(2e)  for  various  Reynolds  numbers: 
□  ,  Re  =  1;  O,  10;  A,  10J;  +,  10s. 

Contour  lines  of  the  components  /  and  g  of  the  axial  velocity  v,/(2«)  as  a  func¬ 
tion  of  radius  r  and  Reynolds  number  Re.  Intervals  are  0.05;  the  tero  level  is 
given  by  the  heavy  line. 

Contour  lines  of  equal  axial  velocity,  v,/(2t)  =  const.,  for  (a)  Re  «=  1;  (b)  10; 

(c)  1 02;  (d)  10s.  Intervals  are  0.01,  0.1,  0.2,  0.2,  respectively.  The  tero  level  is 

given  by  the  heavy  line,  the  velocity  maximum  is  marked  by  +. 

Radial  distribution  of  the  dimensional  velocity  V,  at  r  =  0  for  Re  =  14.9.  The 
symbols  show  the  numerical  solution  to  the  Navier-Stokes  equations  (Vaughn 
1983,  personal  communication).  Parameters:  a  =  60.3  mm,  e/a  =  4.3, 

8  =  20°,  w  =  3000  rpm,  (<  =  500  rpm,  p  =  1400  kg/ms. 

Radial  distribution  of  the  dimensional  velocity  V,  at  t  =  0  for  Re  =  45.7.  The 
symbols  show  the  numerical  solution  to  the  Navier-Stokes  equations  (Vaughn 
1983,  personal  communication).  Parameters:  a  =  60.3  mm,  c/a  =  4.3, 

8  =  20°,  u  =  3000  rpm.  fl  =  500  rpm,  p  =  1400  kg/ro3. 


The  nondimensional  coefficients  rn,,mt  in  eq.  (17)  vs.  ‘he  Reynolds  number  Re 

Comparison  of  the  theoretical  result  for  m,  with:  X  ,  experimental  data  (Miller 
1982)  .O  ,  computational  results  (Vaughn  et  al.  1985).  The  straight  line  shows  the 
asymptotic  law  m,  ~  Re/ 96. 
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Figure  9.  Comparison  of  the  theoretical  reault  for  the  roll  moment  M,  u tu  ■  3000  rpm  vs. 

kinematic  viscosity  v  with:  X,  experimental  data  (Miller  1982)  for 
w  «—  2000  -  4000  rpm;  O,  computational  results  (Vaughn  et  al.  1985)  for 
w  —  3000  rpm.  Parameters:  a  —  60.3  mm,  e/o  «»  4.29,  6  «=  20°,  fl  =  500 
rpm,  p  — ■  1000  kg/m!. 

Figure  10.  Theoretical  results  for  the  roll  moment  M,  vs.  nutation  rate  f)  for  different  spin 
rates:  A,  w  «=  3000;  O,  6000;  0,  8000  rpm.  Parameters:  a  -=  50.4  mm, 
e/a  ■=  4.5,  9  -  20°,  v  —  10*  cSt,  p  —  1000  kg/m*. 

Figure  11.  Theoretical  results  for  the  roll  moment  M,  vs.  spin  rate  u  for  different  kinematic 
viscosities:  □  ,  v  — ■  10s;  O  ,  104,  A  ,  10*  cSt  Parameters:  a  —  50.4  mm, 
e/a  «=  4.5,  9  *»  20°,  0  ■■  625  rpm,  p  —  1400  kg/m*. 

Figure  12.  Experimental  results  for  the  roll  moment  M,  vs.  spin  rate  u  for  different 
kinematic  viscosities:  0  ,  v  •*»  10*;  O  ,  10*,  A  ,  10*  cSt.  Parameters:  a  50.4 
mm,  e/a  •—  4.5,  9  —  20°,  0  «  600  rpm,  p  “  1400  kg/m*. 
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Appendix  B 


His  appendix  give*  listings  of  s  (ample  program,  subroutines,  and  input/output.  Subrou* 
tine  PARAM  converts  dimensional  quantities  as  prompted  into  English  and  SI  units  and  into 
nondimenaional  parameters.  If  the  input  line  starts  with  a  slash  (/),  the  displayed  default  value 
is  used.  If  a  value  of  tero  is  given  for  the  viscosity,  the  program  prompts  for  the  value  of  the 
kinematic  viscosity.  Subroutine  EVALA  calculates  the  functions  /  and  g  in  Appendix  A  .4,  eq. 
(12)  at  r  =*=  k  jK ,  k  «=  0,  ...K  if  K  >0,  as  well  as  J\\)/Rt  and  g\\.)/Re  used  in  Appendix 
A. 4,  eq.  (19)  for  calculating  the  moments.  Subroutine  EVALA  uses  subroutines  BESll  and 
BESIO  in  order  to  calculate  the  ratios  of  Bessel  functions  in  Appendix  A. 4,  eqs.  (12),  (19). 
Additional  information  is  given  by  comments  within  the  programs. 

The  sample  run  provides  results  for  the  run  of  the  Sandiacode  at  Re  14.95  used  in  the 
comparison  Appendix  A.4,  figure  5. 


c  **********  FLUID  FILLED  CYLINDER  *  SAMPLE  PROGRAM  ************* 


C  *  LINEARIZED  EQUATIONS,  Z- INDEPENDENT  FLOW  * 
C  *  WRITTEN  BY  THORWALD  HERBERT,  VP]  &  SU  * 
C  *  * 
C  *  SHOWS  THE  USE  OF  THE  SUBROUTINES  FOR  OBTAINING  A  LIST  OF  * 
C  *  ( 1 )  PARAMETERS  * 
C  *  (2)  VALUES  OF  FIR),  &  <  R ) ,  F‘(1)/RE,  G'(1)/RE  * 
C  #  (3)  DIMENSIONAL  VELOCITIES  AND  MOMENTS  * 


C  *************************************************************** 

c 

COMMON  /FLOW/  REY ,KK ,R< 101 ) ,F( 101 ) , 0 ( 101 ) ,FS1 ,GS1 
COMMON  /ENGL/  RIN,ZIN,OSR,ONR,THD,SG,DENO,VISO, VIKO, VFT 

*  ,MLO,UMOO 

COMMON  /METR/  RMM . ZMM , OS , ON , TH , OK , OY , OZ , OP , DEN , V I S , V I K , VMS 

*  , ML , UMO 
COMMON  /FARM/  AR,TAU,F.FS,RE 
REAL  ML,MLO 
CHARACTER*1  FF 

FF=CHAR( 1Z) 

C 

C***  READ  AND  PRINT  PARAMETERS 
C 

PRINT  2000,  FF 

2000  FORMAT  (A1,’LIN_S  07/25/S5  ***  SAMPLE  PR06RAM ' ) 

CALL  PARAM  (1) 

REY=RE 

C 

C***  SPECIFY  NUMBER  OF  RADIAL  STEPS  KK  (<101 ) 

C 

Kt  =20 
C 

C*ftft  EVALUATE  AND  PRINT  F,G  AT  RADIAL.  POSITIONS 
C 

PRINT  2001,  FF 

2001  FORMAT  (Al) 

H  =MIN(KI  . 100) 

■  'LL  EVALA(l) 

C***  CONVERT  F , G  INTO  DIMENSIONAL  VELOCITY  (M/S) 

C 

PRINT  2002,  FF 

2002  FORMAT  (  Al  r  ,  m ',  7X VELOCITY ,  m/s'/ 

*  15X, ' x, z-p lane  y , z-p I ane ' / ) 

DO  1  1=0, M 

DR  =  RMM*R(  L  +  l )  /  IF? 

DVX=VMS*2 *EPS*F (K+l ) 

DVY  =  VMS*2*EF'S*6  0*1) 

)  PRINT  200?.,  K ,  DR , DVR  ,  DVY 
200'S  FORMAT  <I4,3Fl0.o> 

C 

C***  CONVERT  M  SUB  Z,  M  SUB  Y  INTO  DIMENSIONAL  FORM 
L 

DMZ=-UM0*4*FPS**;*GS1 

DMY=-UM0*2*EPS*(FSl+.25) 

PRINT  2004,  DMZ  , DMY 

2004  FORMAT  (/ 2)'., 'VISCOUS  ROLL  MOW  NT  :  '  ,F  )  1  .6,  '  N  *  rn '  / 

*  2X, 'VISCOUS  PITCH  MOMENT: ' ,FU .6, '  N*m'/i 
PRINT  2001,  FF 

STOP- 

END 
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C  **********  FLUID  FILLED  CYLINDER  ****  PARAMETERS  ************** 


C  *  WRITTEN  BY  THORWALD  HERBERT,  VFI  S  SU  * 
C  *  * 
C  *  SUBROUTINE  FOR  READING  AND  CONVERTING  DIMENSIONAL  VALUES  * 
C  *  INTO  METRIC  UNITS  AND  DIMENSIONLESS  PARAMETERS.  * 
C  *  PRINT  VALUES  IF  I0>0  * 


C  *************************************************************** 

C 

SUBROUTINE  PAR AM  (10) 

COMMON  /ENGL/  RIN,ZIN,OSR,ONR,THD,SG,DENO,VISO,VIKO,VFT 

*  , MLO , UMGO 

COMMON  /METR/  RMM,ZMM,OS,ON,TH,OX,OY,OZ,OP,DEN,VIS,VIK,VMS 

*  , ML , UMO 
COMMON  /PARM/  AR,TAU,EPS,RE 
REAL  ML, MLO 
CHARACTER*4  ID (IB) 

DATA  PI , RADI, RAD2/3. 14159, 57. 2958, 9. 54930/, 1FIRST/1/ 

C 

C***  SET  PHYSICAL  PARAMETERS  TO  DEFAULT  VALUES 

C 

IF  ( IFIRST.EQ.O)  GO  TO  1 

RIN=  2.375 

ZIN=10.375 

0SR=3000 

C'NR=500 

THD=20 

SG= 1 . 4 

DEN=SG*999.b4 
V1S*1 .07E05 
V1K= 1000*VIS/DEN 
C 

C***  READ  ACTUAL  PARAMETERS 

C 

1  PRINT  2000 

2000  FORMAT  (/'PARAM  09/25/84  ***  PARAMETER  CONVERSION '//' ID: 1 ) 
READ  1000,  ID 

1000  FORMAT  (18A4) 

PRINT  2001,  RIN 

2001  FORMAT  (/'RIN  =  radius  (inch)  :  /’,F6.3> 

READ  *,  RIN 

PRINT  2002,  ZIN 

2002  FORMAT  (/'ZIN  =  half-length  (inch)  :  /',F6.3> 

READ  *,  ZIN 

PRINT  2003,  OSR 

2003  FORMAT  (/'OSR  =  spin  rate  (rpm)  :  /’,F6.0) 

READ  *,  OSR 

PRINT  2004,  ONR 

2004  FORMAT  (/'ONR  =  nutation  rate  (rpm):  /’,F6.0) 

READ  #,  ONR 

PRINT  2005,  THD 

2005  FORMAT  (/'THD  =  nutation  angle  (deg):  /’,F6.1) 

READ  *,  THD 

PRINT  2006,  S6 

2006  FORMAT  ( / ' SG  -  specific  gravity  (-):  /',Ft .3) 

READ  *,  SG 

DEN=SG**99.94 
PRINT  2007,  VIS 

2007  FORMAT  (/’VIS  =  v  ,  scc-s  i  ty/zer  o  (cp):  /',F8.0) 

READ  *,  VIS 

VII  -VJ3* 1000/DEN 
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IF  (VIS.GT.O)  60  TO  2 
PRINT  2008 ,  VIK 

2008  FORMAT  (/'VIK  =  kin.  viscosity  (cs):  /',F8.0) 

READ  *,  VIK 
VIS=VIK*DEN/ 1000 
C 

C***  CONVERT  PHYSICAL  PARAMETERS 
C 

2  RFT=RIN/12 

RMM=R  IN*25 . 4 
ZFT=Z  IN/ 12 
ZMM=ZIN#25.4 
AR=7IN/RIN 
OS=OSR/RAD2 
ON=ONR/RAD2 
TAU*ONR/OSR 
TN=THD/RAD1 
0Z=  ON*CO£(TH) 

0Y=  0N*3IN(TH) 

OX=-OY 
OP = OS +07 
EPS=uY/OS 
VFT=RFT*OS 
VMS=VFT#0.3O46 
VEL=VFT-*2*EF'S 
DEN0=1.94*SG 
VlS0=VIS*2.089E-05 
VIk.O=VISO/DENO 
ML*2*PI«RMM**3*AR*DEN*lE-9 
MLO= 2*P I *RFT**3*AR*DEN0 
RE=VFT*RFT /VIKO 
UMO=ML* ( RMM*OS ) **Z* IE -6 
UMOO  =  MLO* ( RFT  *0S ) **2 
IF  ( IO.LE.O )  RETURN 
C 

C***  OUTPUT  SECTION:  LIST  OF  PARAMETERS 
C 

PRINT  2010 

2010  FORMAT  ( / 18 <  ' - '  )/) 

PRINT  2011,ID,RIN,RMM,ZIN,ZMM,AR 

2011  FORMAT  </' ID:  ' ,18A4/ 

*  / 10X, ' A  = 1 ,F10.3, '  IN ' ,F16 .3, ‘  MM' 

*  /10X,'C  = ' ,F10.3, '  IN ' ,F16 . 3, 1  MM',8X,'C/A  =',F10.4> 
PRINT  2012, OSR , OS , ONR , ON , T AU 

2012  FORMAT  ( /9X , 1 03  »',F10.1,'  RPM',F15.2,'  /S' 

*  /9X,'0N  a ' ,F10. 1 , ‘  RPM ' , F 15.2, '  /S',8X,'TAU  =',F10.4) 
PRINT  2013, VFT, VMS. THD,TH, EPS 

2013  FORMAT  ( /8X , ' VEL  ■',FJ0.3,*  FPS',F15.3,'  M/S' //6X, 'THETA  *',F10 

*  '  DEG', FI?. 4,'  RAD' ,7X, 'EPS  =  ' ,F10.4) 

PRINT  20 1 4 , DENO , DEN , S6 , MLO , ML 

2014  FCRMAT  </8X,'DEN  =',0PF10.3,'  SLUG/FT3 ' ,F 1 0 . 1 , '  KG/M3' ,6X, 

*  'S6  =' .F10.3/9X, 'ML  =',0PF10.4,'  SLUG ' ,F 14 . 3 , '  KG’) 

PRINT  2015, VISO, VIS, VIKO, VIK, RE 

2015  FORMAT  ( /5X , 'VIS  =',1PE10.3,'  LDS/FT2 ' ,E1 1 .3, '  CF"  ,/ 

*  /7X, ' V1SK  = ' , 1PE10.3, '  FT2/S' ,E13.3, '  CS',8X,'  RE  =',E10.3) 
PRINT  20l6,UMOO,UMO 

2016  FORMAT  </6X,'U-M0M  = ' , 1 PE  10.3, '  FT-L BS ’ ,E 1 2 .3, '  N*M ' ) 

PRINT  2010 

RETURN 

END 
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C  **********  FLUID  FILLED  CYLINDFR  *  F,G  AND  DERIVATIVES  ******** 


C  *  WRITTEN  BY  THORWALD  HERBERT,  VPI  8,  SU  * 
C  *  * 
C  #  SUBROUTINE  FOR  CALCULATING!  * 
C  *  F  AND  G  AT  KK  EQUIDISTANT  POINTS  IF  KK>0,  * 
C  *  FS1=  F ' < 1 > /RE,  GS1=  G'(1)/RE  * 
C  *  THE  VALUES  ARE  PRINTED  IF  I0>0  * 


C  tt************************************************************** 

C 

SUBROUTINE  EVALA(IO) 

COMMON  /FLOW/  RE ,KK,RR( 101 ) ,FF( 101 ) ,G6< 101 ) ,FS1 ,GS1 
COMPLEX  Z,ZO,S,SO 
DIMENSION  SS(2) 

EQUIVALENCE  (S,SS(1)) 

C 

IF  ( IO.GT.O)  PRINT  2000 

2000  FORMAT  (/'EVALA  09/25/94  ***  FUNCTIONS  F,G  AND  DERIVATIVES') 
C 

C  CALCULATE  F  AND  G 
X0*S0RT(RE/2) 

CALL  BESI1 (XO,ZO,SO) 

IF  (KK.LE.O)  GO  TO  2 
IF  (IO.GT.O)  PRINT  2001 

2001  FORMAT  </3X, ' k ' ,5X, ' r ' ,9X, ' f ' ,9X, ' 9 ' / ) 

K=0 

RR( 1  )=0 
FF  ( 1  )  *0 
GG< l >=0 

IF  (IO.GT.O)  PRINT  2002,  K,RR<K+1> ,FF(K+1 > ,GG(K+1 > 

2002  FORMAT  (I4,3F10.6> 

DR=1 ./KK 

C 

DO  3  K= 1 ,KK 
R=K*DR 
RR(K+1 )=R 

CALL  BESI 1 (R*XO,Z,S) 

S*S*EXP ( Z-ZO ) /SO 
FF(K+1 )=  -SS(2) 

G6(K+1 ) =R-SS( 1 > 

IF  (IO.GT.O)  PRINT  2002,  K,RR(K+1> ,FF(K+1 > ,GG(K+1 ) 

CONTINUE 

EVALUATE  F1 ,G‘  AT  R=1 

CALL  BE3I0(X0,Z,S> 

S='2-S/S0)/RE 
FS1-SS(2) 

GS1=SS( 1 ) 

IF  (IO.GT.O)  PRINT  2003,  FSI,GS1 

2003  FORMAT  ( /5X,9Hf ’ ( 1 ) /Re= ,F10.6/5X,9Hg ' ( 1 ) /Re= ,F10 . 6/ ) 

RETURN 

END 
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c  **********  FLUID  FILLED  CYLINDER  *  BESSEL  FUNCTIONS  *********** 


C  *  WRITTEN  BY  TKCRWALD  HERBERT,  VP  I  «  SU  * 
C  *  * 
C  *  SUBROUTINE  FOR  CALCULATING  THE  MODIFIED  BESSEL  FUNCTION  * 
c  *  ii(2)/E:;p(2) ,  z=>;*u+i>  * 

C  *  * 


C  ******* ******************************************************** 

C 

SUBROUTINE  EESI1(X,Z,S) 

COMPLEX  Z,2Z,S 
DATA  FI  /  '■■ .  1 4 15P265/ 

Z  =  CMPLX(>;,K) 

S=1 

IF  Ol.GE. 10. )  GO  TO  1 

C 

C  ASCENDING  SERIES 
C 

ZZ=Z*I/4 
DO  Z  L=  L ,25 
K=26-L 

2  S=S*ZZ/(K*<K+1 > )+l 
SsS*Z*EXP ( —  Z ) /Z 
RETURN 

C 

C  ASYMPTOTIC  SERIES 
C 

1  ZZ=3*Z 

DO  3  L= 1 ,  IS 
I  =lv-L 

3  S=S*< (2*K-l)**2-4)/(K*ZZ)+l 
S=S/SQRT(2*PI*Z) 

RETURN 

END 
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C  ***#*♦**#♦  FLUID  FILLED  CYLINDER  *  BESSEL  FUNCTIONS  *********** 


C  *  WRITTEN  BY  TNORWALD  HERBERT,  VP1  Si  SU  * 
C  *  * 
C  #  SUBROUTINE  FOR  CALCULATING  THE  MODIFIED  BESSEL  FUNCTION  * 
C  *  Z*JO(Z)/EXP(Z) ,  Z=X*<1  +  I>  * 

C  *  * 


c  ***#**»«■**#****■*»*##**#*******#****■#«•#***■»*»*****#*#*##******«* 

c 

SUBROUTINE  BESIO(X,Z,S) 

COMPLEX  Z,ZZ,S 
DATA  PI/3. 14159Z65/ 

Z=CMPLX(X,X) 

S=1 

IF  (X.GE.IO.)  GO  TO  1 
C 

C.  ASCENDING  SERIES 
C 

ZZ=Z*Z/4 
DO  2  L=l,?5 
K=26-L 

2  S=S*ZZ/(K*K)+1 
S=S*Z*EXP(-Z) 

RETURN 

ASYMPTOTIC  SERIES 

1  ZZ=8*Z 

DO  3  L=  1 , 18 
K=1R-L 

3  S=S*(2*K-1 )**2/(K#ZZ)+l 
S-S*S0RT(Z/(2»PI)> 

RETURN 

END 
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LlN_S  07/25/85  #**  SAMPLE  PROGRAM 
PARAM  09/25/84  ***  PARAMETER  CONVERSION 
ID: 

Sandia  run  at  Re=14.95 

RIN  =  radius  (inch)  :  /  2.375 

/ 

ZIN  =  half-length  (inch)  :  /10.375 

/ 

OSR  =  spin  rate  (rpm)  :  /  3000. 

/ 

ONR  =  nutation  rate  (rpm):  /  500. 

/ 

TiiD  =  nutation  angle  (deg):  /  20.0 

/ 

SG  =  specific  gravity  (-):  /  1.400 

/ 

VIS  =  v i scos i ty/zero  (cp):  /  107000. 

/ 


ID:  Sandia  i 

run  at  Re=14.95 

A 

=  2.375 

IN 

60.325  MM 

C 

=  10.375 

IN 

263.525  MM 

C/A  = 

4 . 3684 

OS 

=  3000.0 

RPM 

314.16  /S 

ON 

=  500.0 

RPM 

52.36  /S 

TAU  = 

0.1667 

VEL 

=  62.177 

FPS 

18.952  M/S 

THETA 

=  20.00 

DEG 

0.3491  RAD 

EPS  = 

0.0570 

DEN 

=  2.716 

SLUG/FT3 

1399.8  KG/ M3 

SG  = 

1.400 

ML 

=  0.5779 

SLUG 

8.434  KG 

VIS 

=  2.235E+00 

LPS/FT2 

1.070F+05  CP 

VISE 

=  8.230E-01 

FT2/S 

7.644E+04  CS 

RE  = 

l . 495E+01 

U-MOM  =  2.234E+03  PT-LBS  3.029E+03  N#M 
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EVALA  09/25/84  *>*  FUNCTIONS  F,G  AND  DERIVATIVES 


k  r  f  3 


0 

0.000000 

0.000000 

0.000000 

1 

0.050000 

0.032883 

0.051961 

2 

0. 100000 

0.065817 

0.103001 

3 

0.150000 

0 . 098828 

0.152194 

4 

0 . 200000 

0.131880 

0. 19§616 

5 

0.250000 

0.164852 

0.241340 

6 

0.300000 

0.197510 

0.279443 

7 

0 . 350000 

0.229474 

0.312012 

8 

0 . 400000 

0.260191 

0.338160 

9 

0.450000 

0.288910 

0.357039 

10 

0.500000 

0.314655 

0.367869 

11 

0.550000 

0.336197 

0.369971 

12 

0. 600000 

0.352040 

0.362801 

13 

0.650000 

0.360400 

0.346007 

14 

0.700000 

0.359201 

0.319480 

15 

0.750000 

0.346067 

0.283427 

16 

0.800000 

0.318339 

0.238448 

17 

0.850000 

0.273092 

0.185618 

18 

0.900000 

0.207178 

0.126587 

19 

0 . 950000 

0.117281 

0.063679 

20 

1.000000 

0.000000 

0.000000 

f 1 ( 1  )/Re  =  -0.176700 
3 1 ( 1)/Re=  -0.034015 
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VELOCITY,  m/5 
x , z-p I ane  y , z - p I ane 


0 

0.000000 

0 . 000000 

0.000000 

1 

0 . 0030 1 6 

0.071047 

0. 1 1226S 

tL 

0 . 006033 

0.142206 

0.222545 

3 

0 . 00c049 

0 . 213529 

0.323333 

4 

0.012065 

0.284941 

0.429133 

5 

0.015081 

0.356182 

0.521444 

6 

0.018098 

0.426743 

0.603769 

7 

0.021114 

0 . 4°5804 

0 .6741 39 

8 

0.024130 

0.562172 

0.730633 

Q 

0.027146 

0.624224 

0.771423 

10 

0.030163 

0.679348 

0.7°4324 

11 

0.03317? 

0 . 726392 

0. 7993o4 

12 

0 .036195 

0.760622 

0.783373 

1? 

0.039211 

0.773686 

0.747537 

14 

0.042227 

0 . 776094 

0.690272 

15 

0.045244 

0.747718 

0.612377 

16 

0.048260 

0.68780S 

0.515194 

17 

0.051276 

0.590047 

0.401049 

13 

0.054293 

0.447632 

0.273506 

19 

0.057309 

0.253399 

0.137586 

20 

0.060325 

0.000000 

0.000000 

VIS 

COUS  ROLL 

MOMENT  J 

3.307994 

N*m 

VISCOUS  FITCH 

MOMENTS  - 

•25.315350 

N*m 
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